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ABSTRACT 

Studying how nuclear star clusters (NSCs) form and how they are related to the growth of the 
central massive black holes (MBHs) and their host galaxies is fundamental for our understanding of 
the evolution of galaxies and the processes that have shaped their central structures. We present the 
results of a semi-analytical galaxy formation model that follows the evolution of dark matter halos 
along merger trees, as well as that of the baryonic components. This model allows us to study the 
evolution of NSCs in a cosmological context, by taking into account the growth of NSCs due to both 
dynamical friction-driven migration of stellar clusters and star formation triggered by infalling gas, 
while also accounting for dynamical heating from (binary) MBHs. We find that in-situ star formation 
contributes a significant fraction (up to ~ 80% ) of the total mass of NSCs in our model. Both 
NSC growth through in-situ star formation and through star cluster migration are found to generate 
NSC " host galaxy scaling correlations that are shallower than the same correlations for MBHs. We 
explore the role of galaxy mergers on the evolution of NSCs, and show that observational data on 
NSC - host galaxy scaling relations provide evidence of partial erosion of NSCs by MBH binaries 
in luminous galaxies. We show that this observational feature is reproduced by our models, and we 
make predictions about the NSC and MBH occupation fraction in galaxies. We conclude by discussing 
several implications for theories of NSC formation. 

Subject headings: galaxies: Milky Way Galaxy- Nuclear Clusters - stellar dynamics 


1. INTRODUCTION 


Over the past two decades, high-resolution observa¬ 
tions with the Hubble space telescope have shown that 
massive stellar clusters reside at the photometric and 
dynamical centers of most interme diate and low lumi¬ 
nosity galaxies of all Hubble types (iCarollo et al.l 119981: 
Matthews et al.l [19991: iBoker et al.l l2002t Balcells et al l 


20031 : iGraham fc Guzmanl 120031 : iCote et al.ll2006D . With 


sizes in the range 2 to 10 pc and masses in the range 
10® to 10® Mq, these nuclear clusters (NCs) have cen¬ 
tral densities up to ^ 10®“’’ Mq pc“®, making them 
the densest stellar systein s observed fe.g.. iPhillips et aP 
rr^lWalcher et a1.ll200,^ . 

NSCs are observed to be larger and brighter, and to 
follow different structural scaling relations than globular 
clusters. The NSC half light radii scale with their total 
mass roughly as ~ while globular clusters have 

rh ~ 3 pc irre spective of th eir luminosity, albeit with a 
large scatter (jHarris 1 [T^ . However, the NSC mass 
distribution overlaps with that of globular clusters at its 
low mass end, suggesting a possible co nnection between 
the two types of s tellar clusters fe.g.. ICote et al.l 120061 : 
[Turner et al.ll201^ . 

NSCs have been studied in detail for large sam¬ 
ples of galaxies in different environments. Ob¬ 
servational work has revealed that the nucleation 
fraction of galaxies is similar for galaxies belong¬ 
ing to very different environments, including the 
Virgo, Fornax, a nd Coma clusters as well as galax- 
ies in the field (|C6te et al.l 120061 : iTnrner et al.l l20l'^ 


Iden Brok et al.l 12014 iBaldassare et al.l 12014) . On aver¬ 
age 80% of all galaxies of all Hubble types with mag¬ 
nitude in the range —20 < Mb —12 contain a well 
defined central cluster. Hence NSCs appear to be ubiq¬ 
uitous structures in galaxies although they tend to be 
missing in galaxies brighter than magnitude Mb ^ —20 
and fainter than Mb ~ —12 (e.g., Ivan den Bergh 1119861 : 
ICote et ^120061) . The real nucleated fraction could be 
however larger than the value constrained from obser¬ 
vations given that the morphological complexity and 
high surface brightness often characterizing the center of 
galaxies can make NSCs difficult to identify, especially 
in massive spheroids. 

High resolution spectroscopic surveys have been used 
to place const raints on NSC age s and star forma¬ 
tion histories (iWalcher et al.l l2006t iRossa et al.l 120061 : 
iBender et al.ll2005f) . The common finding emerging from 
these studies is that NSCs are characterized by a complex 
star formation history with a mixture of morphological 
components and different stellar populations spanning a 
wide range of characteristic ages from 10 Myr to 10 Gyr. 
Cluster ages and masses are also found to depend on 
the host galaxy Hubble type, with NSCs in early-type 
spirals being older and more massive than those of late- 
type spirals. The luminosity-weighted ages of NSCs are 
typically smaller than the ages of host galactic disks, in¬ 
dicating that the growth of the nuclei is a continuous and 
ongoing process occurri ng during and afte r most of the 
host galaxy was formed (|Rossa et al.ll2006ll . 

Age estimates for most NSCs remain however suscep¬ 
tible to large uncertainties due to the fact that the light 
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is often dominated by young stellar populations while at 
the same time the mass is dominated by the old stars. 
The inability to infer spectroscopically the age of old 
cluster components can introduce a bias toward younger 
ages. The bulk of the stellar population is typically found 
in an older and spatially more extended stellar compo¬ 
nent (jCarson et al.ll201^ . This appears to be the case 
for the Milky Way in which more than > 80 % of the 
stars inside the inner parsec formed more than 5 Gyr ago, 
while luminous massive stars dominate the central clus¬ 
ter light within ^ 0.5nc jiB ender et al.ll2005 l: iPfuhl et al.l 
120111 Yusef-Zadeh et al1l2012t iLauer et al.ll2012tl . 

NSC old an d young ste ll ar pop ulations also differ mor¬ 
phologically. iSeth et all (I2006D showed that the three 
edge-on late-type galaxies 1C 5052, NGC 4206, and NGC 
4244 have nuclei that are strongly flattened along the 
plane of their host galaxy disks. Such NSCs show ev¬ 
idence for young stellar components distributed onto a 
disk-like or ring structure superimposed on an older more 
spherical component. In the edge-on spiral galaxy NGC 
4244, old and you ng stellar compon en ts both show clear 
signs of rotation (|Seth et al.ll2008al) . ISeth et all 1)200611 
suggested that the presence of such multiple rotating 
morphological components and the presence of a flat¬ 
tened young cluster aligned with the major axis of the 
host galaxy point toward an in situ build up of nuclear 
clusters, whereby stars form episodically in compact nu¬ 
clear disks and then lose angular momentum or heat 
vertically to form an older spheroidal structure. It has 
been noted however that both rotation and the pres¬ 
ence of young structural components does not exclude 
other formation mechanisms such as for example episodic 
accretion of young star clusters in th e central part o f 
the galaxy due to dynamical friction (|Antoninil l2014f l. 
In the case of NGC 4244 for example, it is likely that 
accreation of star clusters has contributed at least in 
part to the growth of its NSC ()Hartmann et al.l 120111 : 
iDe Lorenzi et alll2013f) . 

Since the early studies, it has been realized that NSC 
masses correlate fairly well with galaxy properties such as 
bulg e velocity dispersion , bulge and galaxy total luminos¬ 
ity (iBalcells et all 120031 : iBoker et al.l 12004 iRossa et al.l 
1200611 . The existence of such correlations suggests that 
the formation of NSCs is intimately connecte d to the for¬ 
matio n and evolution of the host galaxy fe.g.. lLeigh et al.l 
1201511 . Given that NSCs and MBHs are found to coexist 
in some galaxies, that in these systems they have com¬ 
parable masses, and that both follow tight correlations 
with galaxy properties, it is also natural to ask whether 
NSCs and MBHs are somehow connected to each other. 

NSCs and MBHs are known to co-exist in galax¬ 
ies with masses ^ 10^° Mc^ (jSeth et al.l l2008al : 

I Gonzalez Delgado et al.l 1200811 : galaxies with masses 
lower than this value show clear evidence for nu- 
cleation but little evidence for a MBH. Conversely 
galaxies with masses above ~ 10^^ Mq are dom¬ 
inated by MBHs but show no evidence for nucle- 
ation l e.g..lWeh ner fc Harrisll2006HNeumaver fc Walcherl 
I 2 OI 2 II . iFerrarese et al.l ( 2006all foTind that the nuclei and 
MBHs obey a similar scaling relation linking their mass 
to the vi r ial m ass of the host galaxy. More recently, 
iCrahaml ()2012ll obtained NSC - host galaxy correla¬ 
tions using a sample in which he excluded what were 
likely to be nuclear stellar discs from the galaxy sam¬ 


ple of Ferrarese et ah, while including an additional 13 
NSCs in galaxies with velocity dispersions out to about 
200 km/s. Graham found that NSC and MBH masses 
follow different scaling correlations with host galaxy 
properties. Other authors confirmed this result, show¬ 
ing that the mass of NSCs scales with the host-galaxy 
spheroid’s velocity dispersion as Mnsc while the 

mas s of MBHs follows the much steeper relat i on Mmrh 
( iGraham II 2 OI 2 I: lErwin fc: Gadottil 120121: iLeigh et al.l 


2014 Scott fc C rahamI 120131 : iKormendv fc Hd 120131 : 
den Brok et al MU" Although this might suggest that 


MBH s and NSCs did not form from the same mecha¬ 
nism ()Antoninill2013ll . whether the formation of MBHs 
is connected to the evolution of NSCs and whether the 
two types of central objects grow together or in compe¬ 
tition from the same physical process remains unclear. 

The observational findings mentioned above provided 
motivation for theoretical work aimed to understand how 
NSCs form and how their evolution is linked to the evo¬ 
lution of their host galaxy. Two models have been sug¬ 
gested for the formation of NSCs: (i) the cluster infall 
scenario, in which stellar clusters are driven to the galac¬ 
tic nucl eus by dynamical fric tion, merge, and build up 
a NSC (jTremaine et al.l[l975l l: (ii) The nuclear star for¬ 
mation scenario, i n which gas falls into the nucleus and 
forms stars (e.g., iMilosavlievid 12004 iMcLauehlin et al] 
l2006f) . Due to the inherent complexity of gas dynamics 
in star formation, which makes the second of these two 
processes difficult to model, theoretical work has been 
mo stly directed toward th e cluster infall scenario (but 
see lAharon fc Peretsll2015ll . However, both dissipative 
and dissipationless processes are likely to play an impor¬ 
tant role in NSC formation (lAntonini et al.l[2015ll . 

Theoretical studies have employed two different 
methodologies: N-body simulations and semi-analytical 
modeling. 7V-body techniques are used to simulate 
the last stage of inspiral and merger of stellar clus¬ 
ters in the inner region of galaxies. These models 
have shown that a star cluster merger scenario can 
explain without obvious difficulties the observed 
properties of NSCs, including the ir density a n d ve¬ 
locity dispersion profiles (e .g., Bekki et aO 12004 
Caouzzo-Dolcetta fc Miocchil 12008 : lAntonini et ah 


2012 : lAntoninil |201'5 iPerets fc Mastrobuono-Battisti 


2014f) . While these studies make reliable predictions 


about the aspect of NSCs, they suffer from the fact 
that the adopted initial conditions are often not well 
motivated, hence the need to rely on semi-analytical 
models. These are used to make predictions for the 
appearance of star clusters, of known mass and radius, in 
the center of galaxies, and to compute the accumulated 
mass at the center through dynamical friction migration 
over a Hubble time. In semi-analytical models, the rate 
at which the NSC grows by accreting young clusters 
is estimated from empirical cluster formation rates, 
dynamical friction timesc al es an d dissolution tim es (e.g. , 
A^arwaJ_^^Iilosavljeffi ^ I 2 OIII: iGnedin et al.l 12014 


Arca-Sedda fc Caouzzo-Dolcettal 12014) . Studies based 


on semi-analytical approaches have demonstrated that 
the NSC- host galaxy property scaling relations and 
their half-mass radius - mass relation ar e both consisten t 
with formation by star cluster accretion (|Antoninill2013fl . 

All previous calculations assumed NSC formation to 
take place in isolated galactic spheroids, thus neglect- 









































































































formation of galactic nuclei 


3 


ing the role of galaxy evolution, mergers, and the role 
of in-situ star formation. Also, these former idealized 
attempts could not explore the details of the interplay 
between MBH and NSC evolution. In this paper, we 
present a semi-analytical galaxy-formation model that 
allows us to shed light on exactly these points, i.e. it al¬ 
lows us to assess the role of galaxy mergers, MBH merg¬ 
ers and nuclear star formation on the growth of NSCs. 
We follow the formation and evolution of galaxies, MBHs 
and NSCs along cosmic history, including the growth of 
NSCs due to both central migration of stellar clusters 
and in-situ star formation, while also accounting for dy¬ 
namical heating from (binary) MBHs. 

The paper is organized as follows. In Section 2 we 
introduce and discuss the numerical methods employed 
in our study. Section 3 describes the sample data to 
which our numerical results are compared. In Sections 4 
and 5 we describe the main results of our calculations 
and discuss some of their implications in Section 6 . We 
summarize in Section 7. 

2. SEMI-ANALYTICAL MODELS 

We study the formation of NSCs and their co-evolution 
with MBHs along the merger history of their host galax¬ 
ies by semianalytical techniques. 

First, we consider a purely dissipationless formation 
model in which stellar clusters form and migrate to the 
center of a galaxy through dynamical friction. We also 
consider a second semi-analytical model, which unlike the 
first follows the formation and merger history of galaxies 
(both in their Dark-Matter and baryonic components). 
This allows us to study the formation, evolution and dis¬ 
ruption of NSCs in a self-consistent way, accounting not 
only for the dynamical friction-driven inspiral of stellar 
clusters to the nucleus, but also for star formation in 
the nuclear gas, and for the effect of galaxy and MBH 
mergers. 

We begin in this section by describing our methodol¬ 
ogy, highlighting in particular the simplifying approxi¬ 
mations that are made in our models. 

2 . 1 . Cluster inspiral model (CUN) 

In our semi-analytical model CUN a galactic nucleus 
forms through the orbital decay (via dynamical friction) 
and merger of star clusters in the central region of a 
galaxy. This code adopts a rather idealized model of 
an isolated galactic spheroid and follows the consecutive 
inspiral of stellar clusters along their orbits until they 
decay into the galactic center or they are disrupted in 
the process. This relatively simple approach has the ad¬ 
vantage that the details of the tidal interaction of the 
clusters with the background galaxy and central MBH 
are easily included, and in a relatively precise manner. 
On the other hand, as previously stated, this model has 
the important limitation that galaxy evolution as well 
as dissipative processes that can lead to star formation 
episodes in the galactic center are neglected. Another ba¬ 
sic assumption made in the CUN model is that MBHs are 
already in place before NSCs grow around them; more¬ 
over we assume that all galaxies contain initially a MBH. 

The code is essentially the same as that used by 
lAntoninil (|2013f) . and we direct the reader to that pa¬ 
per for a more detailed description. Briefly, we represent 
the galaxy spheroid by a simple power-law density model; 


p{r) = po (r/ro) , where po = {3 - 7 )Msph/ 47 rri^. This 
latter expression a ssumes t h at the density of the galac¬ 
tic bulge follows a iDehnenI (|1993D profile, r ~ r~^ , at 
r ^ tq. Given a mass for the central MBH, Mmbh, we 
compute the galaxy density profile slope, scale radius, 
spheroidal mass and galaxy velocity dispersion by using 
the corresponding scaling relations tha t link these prop¬ 
erties to Mmhh (iGiiltekin et al.l 120091 iGrahamI l2012bl 
iGraham fc Driver] 1200711 . The galaxy effective radius, 
ReS, was derive d from the s ize-ma ss relation given by 
equation (32) of IShen et al.l (|2003ll . T he ratio Resjrp 
follows from equations (4) and (17) of iDehn^ (Il993f) . 
which give Res/ro = (1.8, 1.5, 1) for 7 = (1, 1.5, 2). 

We assign the stellar cluster masses using the clus¬ 
ter initial mass function (CIMF), dn/dnigc oc. rn~^ 
(|Bik et al. I 120031 Ide Griis et al. I I2003H . and limiting 
mass values of mmin = IO^Mq, rumax = 1 O^M 0 . The 
stellar clusters are assumed to form continuously over 
the age of the galaxy. 

As commonly done, we assume that the clusters have 
initially the same distribution as the stars in the galaxy, 
and that initially a fixed fraction /gc = 0.05 of stars 
form in strongly bound stellar clusters. This latter is ap¬ 
proximately consistent with the typical cluster formation 
efficiency for the Milky Way found by iKruiissenI l) 2012 [) 
and it is similar to the value adopted in previous stud¬ 
ies (|Gnedin et al.l 1201^ . Gnedin et al. showed that by 
adopting a power-law mass function dn/drrigc oc m~^ 
and a fraction /gc = 0.04 of clustered star formation, one 
can reproduce the observed density profile of clusters in 
M87. 

We assume that the central properties of a stellar clus¬ 
ter remain unchanged during inspiral and that rt > rx, 
where tk is the cluster core ra dius and rt i s the cluster 
tidal (limiting) radius given by (|King 1119^ : 


n 


<JK 





( 1 ) 


with (j) the sum of galactic, MBH and NSG potentials, 
and a a “form factor” that depends on the density dis¬ 
tribution within the cluster. The mass of a King model 
is then related to its tidal radius and velocity dispersion 
via the expression 


Gmgc ~ . (2) 

We note that although at first order the central prop¬ 
erties of the clusters are not expected to change during 
inspiral, the stars in a NSG formed from cluster infalls 
should arguably end up with the same root-mean-square 
velocity as the galaxy host , which appears to be consis¬ 
tent with observations Isee lLeigh et al.ir2ni5h . 

Given a cluster of central velocity dispersion aK, the 
time evolution of its orbital radius is 

1 

3-7 

(3) 


-(t) = 


3 — 7 
r- 


(3-7^ 


(4 - 7 ) 1/27 47rGpo 


F{-f) In Aal x t 


where is the initial separation from the center, for 
the Goulomb logarithm we take In A = 6, an d F{-f) is 
given in equation (1) of iMerritt et al. I (1200411 . We ap¬ 
proximately account for the dynamical dissolution of the 
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Accretion and mergers Dynamical friction. 



Fig . 1. — Schematic representation of the model of IBaraussel 11201211 (in black; including the improvements described in ISesana et al.l 
Hail)) with the additions described in this paper (in red) to model the formation and evolution of NSCs. 


clusters due to their collisional evolution and interaction 
with the external tidal field of the galaxy by only evolv¬ 
ing clusters with dynamical friction time shorter than 
their dissolution time [equation (fTTI) below]. 

The cluster orbits are evolved using equation ([31) up 
to a maximum integration time of lOGyr, and the NSC 
mass is obtained by summing up all the mass transported 
by the infalling clusters within a radius tnsc- The NSC 
radius, rNSCi is computed fro m the recursive relation 
given by equations (14)-(16) of I Antoninil (I2013D . As we 
evolve the cluster orbits, we compute their tidal disrup¬ 
tion radius due to the galaxy, MBH and a pre-existing 
NSC gravitational fields, adopting a cluster core radius 
tk = Ipc, roughly equal to the median value of the core 
radii list ed in the Ha rris’ compilation of Galactic star 
clusters (|Harris I [l996D . When the clusters reach their 
tidal disruption radius, the remaining core mass is dis¬ 
persed around the radius of disruption. When evaluat¬ 
ing the final NSC density distribution we make the as¬ 
sumption that the stars from the disrupted clusters are 
isotropically distributed around the center. We note that 
in reality the clusters will distribute initially in a disk¬ 


like structure, but this will morph into a more spheri¬ 
cal distrib ution over a fra ction of the nuclear relaxation 
timescale (IAntoninill2014l) . 

A MBH binary formed during a galaxy merger leaves 
an imprint on the galactic nucleus in the form of a 
mass deficit, Mej;bin, a decrease in the mass of the 
nucle us due to ejection of stars on intersecting or¬ 
bits (jBekki fc CrahanTI I2010D . We account for this by 
subtr acting a mass Mej;bin ~ 0.5 x N^Mmbh (|Merrittl 
|200(3D . where is the predicted number of major merg¬ 
ers after the last major accretion event, from the final ac¬ 
cumulated mass inside tnsc • The number of mergers 
is obtained from the gal axy mass using the distribution s 
shown in Figure (2) of iHaehnelt fc Kauffmaniil (I2002D . 
where in order to convert H—luminosity to mass we used 
a constant mass-to-light ratio M/Ly = ^Mq/Lq^ as 
sometimes adop ted in the literatur e for an old stellar 
population (e.g.. iCnedin et ai1l2014D . 

Finally we compare the resulting NSC density profile, 
obtained after 10 Gyr of evolution, to the density profile, 
p(r), of the galaxy, and define a galaxy as nucleated if 
at some radius the NSC stellar density was found to be 
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higher than that of the background galaxy. This latter 
condition is based on the fact that if the NSC densities 
are below the galaxy density at all radii, then it would 
be difficult to observationally distinguish the two com¬ 
ponents. 


2.2. Galaxy-formation model (GxeV) 

Our galax y formatio n mode l (hereafter GxeV) is based 
on that of iBaraussd (20121. which was fu rther im¬ 
proved in iSesana et al.l (j20l4) . More precisely, iBaraussej 
lj2012l l presented a semi-analytical galaxy formation 
model tracking the evolution of baryonic structures 
along Dark Matter merger trees. These trees are pro¬ 
duced with an extended Press-Schechter formalism, mod¬ 
ified in order to reproduce the results of N-body sim - 
ulations of Dark-Matter halos (jParkinson et al.l 1200811 . 
while the baryonic structures include the hot, largely 
unprocessed intergalactic medium, galactic disks and 
spheroids (in both their stellar and cold, chemically- 
enriched interstellar medium components), a low angu¬ 
lar momentum reservoir of cold nu clear gas available for 
accretion onto t he central MBH (jCranato et al.l 12004 
iLani et HI 12014) . and the MBH itself. These baryonic 
components are interconnected by a plethora of interac¬ 
tions, schematically summarized in Figure [TJ Also in¬ 
cluded are environmental and tidal interactions between 
galaxies during mergers, following simple semi-analytica l 
recipes (iTaffoni et ai.ll20 0,4 iBovlan-Kolchin et al.l[2008[l . 
Moreover. ISesana et al.l (j20l4l later improved this 
model by refining the star formation law - adding 
in particular an explicit dependence on the metallic- 
ity (|Krumholz et al.l [200911 . as well as a distinction be¬ 
tween pseudo-bulges forming from bar instabilities and 
classical bulges created by major mergers - and by devis¬ 
ing more realistic prescriptions for the evolution of the 
spins of MBHs with cosmic time. 

For this work, we have further ameliorated the model 
by including the formation and evolution of NSCs (cf. 
Figure [T|). More precisely, we consider two possible for¬ 
mation channels for these objects, namely one where 
NSCs form from the dynamical friction-driven migration 
to the nuclear regions of star clusters created further out 
in the galaxy, and one in which NSCs form in situ as 
a result of star formation in the nuclear regions. Note 
that we assume no high redshift seeds for the NSCs, i.e. 
the NSC mass grows from zero at high redshifts to its 
present-day value through the two aforementioned chan¬ 
nels (this growth being modulated/hindered by the effect 
of galaxy and black-hole mergers, cf. section [2. 2. 3D . 

In the GxeV models described below we define as nu¬ 
cleated any galaxy containing a central cluster more mas¬ 
sive than 10^ AIq, as this value corresponds approxi¬ 
mately to the l ower limit of the obse rved NSC mass dis¬ 
tribution le.g.. lden Brok et al.ll2ni4l . We also define an 
early- (late-)type galaxy as one with bulge to total mass 
ratio B/T > 0.7 (< 0.7). For each of the GxeV models 
that we present, we simulate about 1300 galaxies with 
Dark-Matter masses ranging from IO^^Mq to IO^^Mq. 

2.2.1. Migration channel 

To implement the dynamical friction-driven migra¬ 
tion of star clusters to the central nucleus, we first as¬ 
sume that star clusters are created during star-formation 


events with efficiency fsc, bot h in galacti c bulges and 
disks. In bulges, IBaraussd l)2012ll and ISesana et al.l 
(|20l4) assumed a volumetric star-formation rate 
while in disks they considered a vertically-averaged rate 
of star formation S,^. Therefore, we assume that star 
clusters form in bulges and disks, and that the total 
mass of the star-cluster population (in both the bulge 
and disk), Mgc, changes with rate 


-^formation ~ /gc ^ 


47r / Pi.r'^Ar -\-2Tr i]*rdr I , (4) 


where we choose fp-r = 0.07, whi ch is the typical value in 
the Milky Way (jKruiissenll2012li . We stress that plausi¬ 
ble values /gc = 0.05 - 0.2 would only mildly impact the 
normalization of our results for the Mnsc ~ o' relation, 
and would not affect its slope significantly. In addition 
we have tried a model with a variable /gc set to 0.07, 0.04 
and 0.5 in disk, quiescent and starburst galaxies respec¬ 
tively; these correspond app roximately to the values of 
Figure 4 in iKruiissenI (j2012[ ). The results of this model 
were found to not significantly differ from those of mod¬ 
els with a fixed /gc. Thus, for the sake of simplicity, we 
present in what follows only the results of a model with 
fixed /gc = 0.07. 

We then assume that the star clusters formed in the 
bulge are spatially distributed in the same way as the 
bulge stellar population, i.e. if that population has 
density profile Pbuige(^)i volumetric number density 
of star clusters is assumed to satisfy 




ptuige(^) 


(5) 


where A'gc'^^®® is the total number of star clusters in the 
bulge, and Mb „ige is the total bulge mass in stars. Note 
that following iBarauss i dml, we assume a Hernq uist 
distribution for Pbuige(^); refer to IBaraussd (|2012l ) for 
more details on the choice of the scaling radius for this 
density profile. 

Similarly, we assume that the disk’s star clusters are 
distributed according to the same surface density profile 
as the overall stellar content of the disk, i.e. the superfi¬ 
cial number density r'gc’^''(r’) will satisfy 


P 


disk 

gc 


(r) = 


/ydisk 
^ gc 


-^disk 


( 6 ) 


where again is the total number of star clusters in 

the disk, and is the mass of the s tellar disk. Again, 

for Dj;gj.(r) we follo w IBaraussd (1201211 and adopt an ex¬ 
ponential profile Isee lBaraussd (|2012[1 for more details on 
the calculation of the scaling radius for this distribution). 

As for the mass distribution of the newly formed 
star clusters, we assume a power-law mass function 
dp^j./dmgc oc TTigc”^ (with X = bulge, disk). Combined 
with equations m and ([6]), this gives the distribution 
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functions for the bulge and disk star clusters 


star clusters in isolation through 


, , 1 d^iV^"^se 

.bulge _ ^ _ gc 

drUgcdV 


dm 


gc 


‘'gc 


TT 


disk _ 
gc = 


1 

dmgcdS' 


^^ PscHr) 

dmgc 


(7) 

( 8 ) 


where dl^ and dS” are units of volume and surface. We 
assume that these distribution functions are valid for 
masses Wgc for individual star clusters between mmin = 
IO^Mq and mmax = 10 ® - IO^Mq, and the normalization 
factor 


A = 


^max ^min 


^max ^min 


(9) 


ensures that the integral over all masses and over the 
whole volume (or surface) is 1. Note that in our runs we 
assume mmax = IO^ Mq; using the dissipationless model 
described in Section 12.11 we find that reducing mmax to 
1 O®M 0 has the effect of changing the normalization of 
the scaling relations by a factor ~ 10 , but did not change 
their overall slopes. 

These distribution functions can then be used to cal¬ 
culate approximately how much mass is lost by the pop¬ 
ulation of star clusters due to their collisional evolution 
and interaction with the external tidal field of the galaxy 
(“dynamical dissolution”): 


dissolution 
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( 10 ) 


with 


ttid ~ lOGyr 


mcl 


2 X lO^M© 


P{r) 


( 11 ) 


and 


P{r) = 41.4 


Kir 


kpc J \km/s 


( 12 ) 


Kir b eing the halo’s virial velocity (e.g., iGnedin et al.l 
[Ml. Recent iV-body sim ulations show that a ~ 2/3 
([Gieles fc Baumgardtll200l . which we adopt here. Also, 
note that for simplicity we do not track separately the 
mass in star clusters in the bulge and that in the disk, 
but simply follow the evolution of the total mass in star 
clusters Mgc [cf. equation (|1|)]. In equation (ITUl) [and 
in equations (fOl) and (ITBl) below] we therefore simply 
assume that the mass of star clusters in the disk is ~ 
[■^disk/(117buige + mdisk)]-^gc: and that the mass of star 
clusters in the bulge is « [M^uige/iM^nise + M^isk)Wgc- 
In a similar fashion, we account for the evaporation of 
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where 


17md 


2 X 10® M© 


Gyr. 


(14) 


Note that because the evaporation time-scale does not 
depend on the star cluster’s location (since it describes 
evaporation in isolation), equation (IT^ can be simplified 
to 


-^evaporation ~ ^-^gc 


r-2 

'•gc 


tev (md) 


dmd. (15) 


Note that our models do not include the possibility of 
close encounters with giant molecular clouds, which in 
some cases could significantly reduce the clusters’ life¬ 
time. Also they rely on necessarily simplihed models for 
the galaxy. A more realistic model for the galactic po¬ 
tential might somewhat change the timescales introduced 
above. However, we note that analytical models of NSG 
formation indicate that the final NSG mass is not much 
affected by the assumed slope o f the galactic ba ckground 
density profile (see Section 4 of IAntoninill2013l) . 

The total mass in star clusters also decreases because 
the individual star clusters fall toward the central NSG, 
under the effect of dynamical friction, i.e. 
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(16) 

The dynamical-friction timescale in the bulge is domi¬ 
nated by the interaction with the stellar background, and 
we therefore follow Binney and Tremaine (2008) [equa¬ 
tion ( 8 . 12 )]: 


idf,buige(A«T-d) ~ 15Gyr- 
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10^Mf7 


log A \5kpcy lOOkm/s md 

(17) 

where a is the bulge velocity dispersio n (related to th e 
halo’s virial velocity by tr ~ 0.65Kir ([Ferraresd I 2002 IH 
and we choose log A = 6 . As for the dynamical friction 
in the disk, we account for both the interaction with the 
stellar and gas component, i.e. 


idf,disk(?',md) = idf.disk(AWd) +iIf"disk(Amd) 
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where L is the orbital angular moment of a star clus¬ 
ter of mass md belonging to the disk (and thus moving 
with velocity ~ at a separation r from the galactic 
center), while the angular mo mentum loss due to dy¬ 
nami cal friction is expressed (|Agarwal k, Milosavlievid 
[miD in terms of the surface densities S* and Sgas 
of the disk’s stellar and gaseous components, the ve¬ 
locity dispersion of the stars in the disk ct*, and the 
velocity dispers ion/sound speed of the gas in the disk 
o^gas ~ O.lcr* (iDuttoii fc van d en Boschl p009f) . Note 
that for cr* we follow iKreeel et alJ ( 2005!) arid assume 
cr* K, 0.29t4nax, where Vmax is the maximum rotational 
velocity inside the disk. Also observe that in equation 
Cl we are summing the inverses of the timescales, so 
that the infall rates due to the dynamical friction from 
the gas and that due to the stars get summed in equation 

Cl. 

Not all of the mass of the star clusters falling toward 
the NSC eventually accretes onto it, because star clusters 
get truncated due to the galactic tidal field, and tidally 
disrupted by the central MBH. To account for these ef¬ 
fects we assume that the rate of change of the NSC mass 
due to infall of star clusters is 


= ( 21 ) 

where clearly the major difficulty lies in computing the 
factor 0 < < 1. Since it would be computationally pro¬ 

hibitive to evolve the infall of the individual star clusters 
self-consistently within our semi-analytical galaxy forma¬ 
tion model, we calculate F with the “monolithic” CUN 
model of Section EH in order to derive an easy-to-use 
analytical expression. 

More specifically, approximating a star cluster as a 
King model, its tidal radius at a distance tnsc from the 
center of a galaxy containing a MBH and a NSC at its 
center is: 
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V2 


47rGpo 


-k 


3G(Mcmo) 


'NSC 


J'o / 3-7 

- 1/2 


( 22 ) 


where Mqmo = /Hmbh + Mnsc is the total mass of the 
central massive objects (i.e., MBH plus NSC mass), and 
for the galaxy density profile we adopted the power-law 
model p{r) = po (r/rp)”^. Note that to be consistent 
with the Hernquist profile used in GxeV, we set 7 = 1 in 
equation (l 22 l) . Also, the distance tnsc is to be set to the 
outer radius of the NSC, i.e. the distance from the galaxy 
center below which the star cluster is assumed to have 
become part of the NSC. We assume r-NSC = 5rh, where 
r-h is the NSC half-mass radius. Our detailed prescrip¬ 
tion for Th will be presented in equation (1251) below. For 
a King model, the truncated mass mgc(fT_R-) of the star 
cluster is then related to its tidal radius via equation ©• 
Using the CUN model described in Section 12.11 we 
find that more than 90% of the NSC mass comes from 
star clusters with initial mass > O.lmmax, and that the 
(mass-weighted) mean central velocity dispersion of the 
star clusters contributing to the nucleus growth is {(tk) ~ 
20 km/s. When setting mmax = 10^ Mq, the average 
initial mass of these star clusters is {rrigc,in) ~ 2.5 x 
10 ®Mq. Assuming that all star clusters that decay to 


the center have similar central properties {ctk)-, and that 
these properties do not change during the infall, we have 


Wgc((g/c)) 
(777-gc, in) 


(23) 


where TOgc((o’K)) is given by equations (E)) and (EH)- Fi¬ 
nally, to account for the possibility that a star cluster 
may undergo complete tidal disruption before it decays 
to a radius t-nsCi we compare the tidal radius given by 
equation (IMl) with the core radius rjf « 1 pc of the star 
cluster. If rt < r^, we set = 0. 


2.2.2. Star-formation channel 

As mentioned above, NSCs may also form by star 
formation in the cold gas accumulating in t he galac¬ 
tic cen ter. In the galaxy formation model of iBaraussd 
(1201211 ('and in t h e imp roved version of this model used 
bv iSesana et al.l (12014 11. transfer of cold gas to a low- 
angular momentum nuclear reservoir - available for ac¬ 
cretion onto the central MBH - is assumed to be cor¬ 
related with star-f ormation events in the galactic bulge 
through a relation (Granato et al.ll 20 dl: iLaoi et aI1l20l4 


iHaiman et al.l[2f)f)4) 


M/Xn = AesV’b , (24) 

where ' 0 b is the star formation rate in the gaseous 
bulge, and Ares ~ 10 “^ - 10 “^ is a free parameter, 
which we set to A r^s ~ 6 x 10 “^ in t his paper (as in 
ISesana et al.l (l201-(^l. I n the model of|Baraussi (|2012f l 
and ISesana et al.l (|2014ll . star formation in the bulge is 
associated with major galactic mergers, and (less impor¬ 
tantly) with bar instabilities in the galactic disk. 

Note that the physical mechanism responsible for the 
loss of angular momentum of the cold gas and its 
transfer to this nuclear reservoir may be the radiation 
drag caused by the newly formed stars (IU memur all 2001 l : 
iKawakatu fc Umemural[200^ iKawakatu et al.l 1200311 . or 
the reshuflling/shocks of the gas as a consequence of disk 
instabilities or mergers. More generally, a correlation be¬ 
tween bulge star formation and MBH accretion/growth 
is expected based on the M — a relation for MBHs and on 
the parallelism be tween the quasar luminosity and star 
formation history (iGranato et al.l[2004 ILaoi et al.l[20l4 

IHaiman et al.l[2004b 

The reservoir’s gas is then made available to accrete 
onto the MBH on the viscous timescale. Whenever 
the resulting viscous accretion rate exceeds AEddAfEdd 
(-Mfidd being the Eddington mass accretion rate and 
Asdd a free parameter), we truncate the accretion rate 
to that value. To allow for moderately super-Eddington 
accretion in the case in which MBHs form from light 
pop HI star seeds (of ma ss ~ 200 Mq at 2 ; ~ 15 — 
20) (iMadau fc Reea [200111 . we set AEdd = 2.2 for that 
seed model. This is because some amount of super- 
Eddington accretion is known to be needed to recon¬ 
cile light MBH se eds with the quasar luminosity function 
at high redshift (|Madau et al.ll20l'3 l. We also consider 
two “h eavy” seed m odels (namely that o fiVolonteri et al.l 
(|2008f) and that of iKoushiappas et al.l (|2004lH whereby 
MBHs grow from seeds of ^ IO^Mq at 2 ; ^ 15 — 20, in 
which case we set AEdd = 1 - 

On top of these prescriptions, in this paper we assume 
that the gas in the nuclear reservoir also forms stars. 
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which are assumed to contribute to the NSC. To com¬ 
pute the star-formation rate, we need to choose a size 
for the reservoir. It is natural to assume that this size is 
comparable to the observed size of NSCs. More specifi¬ 
cally, we assume that the reservoir is disk-like, and has a 
(vertically-averaged) exponential density profile, whose 
scale radius we choose such that the reservoir has the 
same half-mass radius Vh as the NSC. The latter is as¬ 
sumed to be 


We can thus summarize the depen dence of fr on the 
metallicity by the exp licit expression (iSesana et al.ll2014 
iKrumholz et al.ir200^ 



( 


fc=< 

1- 

1 + )™) 


, 2%, 



if /c > 2%, 

otherwise 


(27) 


with 

s = ln(l -b0.6x)/(0.04EiZ'), 


X = 0.77{1 + 3.1Z'°-^^^), 


with Mdyn = Afres + M^sc the dynamical mass of the 
nucleus (including the gas in the nuclear region). Note 
that this scaling is inspired by observations of the size of 
NSCs, whic h is found to scale w ith the square root of the 
luminosity (| Turner et al .11201^ . 

Once a density profile for the reservoir is specified, 
we apply a sta r form ation law similar to that used 
bv ISesana et al.l (1201411 for gal actic disks. No te that 
ISesana et al.l ( 20141) improved on lBaraus^ (120121) by con¬ 
sidering different modes of star formation in classical 
bulges forming from major mergers - where star forma¬ 
tion is assumed to take place in “bursts” of duration 
comparable to the bulge dynamical time - as opposed 
to disks and pseudobulges arising from bar-instabilities 
of disks - where star formation is assumed to happen 
in a “quiescent” mode , described by the prescription of 
IKrumholz et al.l (120091) . Since star formation in the cen¬ 
tral region of our Galaxy is known to be weaker by at 
least a factor ten relative to what wo uld be expected 
based on the observed gas densities (jKruiissen et all 
120141) . it seems more appropriate to use the “quiescent” 
star-formatio n mode for ou r nucle ar re servoir El Follow¬ 
i ng th erefore ISesana et al.l (j2014l ) and IKrumholz et all 
(l2^ . we assume that the star formation in the reser¬ 
voir takes place on a timescale tsF and involves only a 
fraction /c of the reservoir’s gas, thus allowing us to write 
an approximate star formation rate in the reservoir as 

= ( 26 ) 

CSF 

Note that this simple prescription does not account for 
the spatial distribution of the star formation activity, e.g. 
the observational fact that younger stellar components 
seem to be more centrall y concentrated than the r est of 
the NSC [see for example iGeorgiev k. Bokeil (|2014D ]. 

The fraction /c of cold gas available for star formation 
is set by the fraction of molecular gas for metallicities 
Z' > 0.01 (in solar units), when star-formation happens 
in molecular clouds. At lower metallicities Z' < 0.01, 
star f ormation takes place in the atomic phase ([Krumholzi 
120121) . In general, /c is expected to decrease as the metal¬ 
licity decrea ses, but recent ob servations of nearby spirals 
and dwarfs (iBigiel et al.||2010D, as well as the Small Mag¬ 
ellanic Cloud (IB olatto e t al.ll20lH) show that it levels off 
at 2% ([Genzel et al.ll2015l) . 

^ Observe that although star formation takes place quiesc ently 
in our model’s nuclear reservoir, we also assume [cf. equation (I24ll 1 
that the reservoir’s feeding is triggered by star-formation events in 
the bulge, i.e. mainly by turbulence-driven bursts of star formation 
following major mergers. 


(5 = 0.0712 (0.1s“^-f 0.675) , 


El = Eres/(M0 pc ^) . 

The timescale fsF is in stead given by (jSesana et al.ll201^ 
IKrumholz et al]r2009D 


tgp = (2.6 Gyr) ^ x 


/ s-0.33 

j , Eres > Sth 


(28) 


with Eth = 85M0/pc^. This expression is obtained by 
assuming that star formation happens in clouds, and 
the two branches appear according to whether the cloud 
density is determined by internal processes, or by the 
external pressure (in galaxies with sufficiently high sur¬ 
face densities). Also, note that the typical cloud mass 
does appear in this expression because it is related to 
the gas surface density by identifying it with the local 
Jeans mass, and by assuming marginal gravitational sta¬ 
bility of the reservoir. For the surface density of the 
reservoir, in both equation and (l2^ . we choose 
Eres ~ Afi.es/(27rr^gg), Fres being the scale radius of the 
reservoir’s exponential surface density profile. This is 
the central value of the density for an exponential sur¬ 
face density profile. While this choice is quite arbi¬ 
trary, we note that at high surface densities the com¬ 
bination fc/tsF lias a weak dependence on Ej-es, namely 
fc/tsF oc E°gg^. Therefore, a different choice of Ej-es by a 
factor 2 (5) only changes as given by equation (1^ 
by 0.1 (0.24) dex. 

We stress that equation (1^ is admittedly a rough pre¬ 
scription for the star formation in the nuclear regions. 
However, it seems to work reasonably well when com¬ 
pared to observations of the star formation rates in the 
central 500 pc of the Milky Way (the “central molecular 
zone”). Indeed, in Figure|2]we compare the predictions of 
our star formation prescription [where in equations (I27|) 
and (1^51) we identify Ej-es with the surface gas density 
represented in the horizontal axis] with measurements of 
the star formation in the central molecular zone at differ¬ 
ent angular scales (and thus different average gas surface 
densities). 

In the light of this reasonable agreement, we chose not 
to include any feedback effects (AGNs or supernovae) 
when considering star formation in the nuclear regions. 
We stress however that both kinds of feedback are in¬ 
cluded when consideri ng the bulge and disk star forma¬ 
tion [cf. Figure [U and iBaraussd (|20il) ] . They therefore 
indirectly affect the star formation in the nuclear region. 
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Fig. 2.— A comparison between our star-formation prescrip¬ 
tion in the nuclear region — denoted by a blue line, the width of 
the line representing the scatter of the predicted star-formation 
rate with metallicity, which we assume in this figure to vary be¬ 
tween 1 and 4 solar units - and the measured star formation in 
the Milky Way ’ s cen tral molecular zone (CMZ) as reported by 
IKruiissen et a U ll20Tl'l : the average star formation rat e in the cen¬ 
tral 1 50 pc of NGC 253, as reported in Table 6 of I Leroy et al.l 
lf20T5l) : the Milky Way’s star formation rate at galactocentric dis¬ 
tance of 1 kpc and that of NGC 69 46 at galactocentric distan ces 
of ~ 300 and 700 pc, as reported bv IKennicutt &: Evan^ l|2012l') in 
their Figure 7; and the measured nuclear star formation in the star- 
forming galaxy SDP 81 at 2 3 l|Dve et al.ll2015l : ISwinbank et all 

1201511 . For comparison, we also show by a dashed orange line the 
prediction of the Kennicutt-Schmidt law IIKennicuttlll998) . 

since the feeding of the reservoir depends on the larger- 
scale galactic evolution via eqnation (IMl) . 

2.2.3. Effect of galaxy mergers and black-hole binaries 

Delays between galaxy and black-hole mergers. — 

When two dark matter halos merge, the two galaxies, 
being more compact than the halos, initially retain their 
identity, and are slowly brought together by dynami¬ 
cal friction. During this evolution, environmental effects 
such as tidal stripping and tidal evaporation become im¬ 
portant and contribute to remove mass from the smaller 
galaxyQ When the galaxies hnally merge, within a few 
dynamical-friction time-scales after the halo merger, the 
MBHs of the two galaxies are expected to form a binary 
system. 

This binary is slowly brought together by a variety 
of effects, i.e. dynamical friction at large separations, 
and at smaller < pc separations by gas interactions (if 
sufficient gas is present in the nuclear regions) or by 
interaction with stars if the nnclear star cluster has a 
sufficiently dense core, or by triaxiality of the central 
potential. Gas interactions (and in particular planet¬ 
like migration) might bring a MBH binary to merger on 
time- scales < 10^ — 10® yr (|Haiman et al.l 120091 iColoil 
120141 , while three-body interacti ons with stars might re¬ 
quire timescales up to ^ Gyr (|Yul 1200211 or more. In 
fact, it is not completely obvious that a MBH binary 
would merge at all (within a Hubble time) in the pres¬ 
ence of three-body interactions with st ars alone. This 
is kn own as the “final-parsec” problem (jBegelman et al.l 
llOSOtl . but recent numerical simulations suggest that in 
triaxial merger remnants (such as those that would be 
expected from a recent galaxy merger) MBH binaries 

^ We refer to IBarauss i Oom for more details about how we 
compute the dynamical-friction timescale and these environmental 
effects. 


migh t coalesce on timescales of a few Gy r (iKhan et al.l 
120111 iVasilievI 12014 iVasiliev et al.l 12014 1201511 . Ro- 
tation in flattened galaxy models has also been sug- 
gested to help drive MBH binarie s to coalescence, cf. 
iHollev-Bockelmann fc KhanI (l2015tl . Also, if a MBH bi¬ 
nary stalls at ~ pc separation, a later galaxy merger 
may add another MBH to the system, forming a MBH 
triple. At least for large ^ 10®Mq MBH masses, triple 
MBH interactions may trigger binary mergers (and pos- 
sibl y eiection of the smalle st MBH) on timescales < 10® 
yr (|Hoffman fc LoebI 1200711 . Because of the uncertainty 
about the mechanisms described above, we adopt here 
a set of minimal simplified prescriptions to estimate the 
“delay” between galaxy and MBH mergers. These pre¬ 
scriptions are presented in the following. 

In gas-rich environments, defined by the criterion 
Mres > A1 mbh,i + A 1 mbh, 2 , we assume that the de¬ 
lay is simply given by the viscous time of the nuclear 
gas. At a distance r from the central MBH, the vis¬ 
cous time is simply ti, ~ jv^ v being the kinematic 
viscosity. The latter can be approximately expressed in 
terms of the gas velocity v and the Reynolds numbe r 
TZ as V rv/TZ ([Frank et al.l 120021 iDuschl et al.ll200^ . 
Since the kinematic viscosity is huge for accretion disks, 
it is natural to assume that TZ is the critical Reynolds 
number 7^c - 10^ - 10® marking the onset of turbu¬ 
lence. Indeed, laboratory experiments find that once 
this critical Reynolds number is reached , the viscosity 
increases significantly ([Frank et al.l l2002t iDuschl et al.l 
120001 1. Assuming now that r ~ G(MMBH,i+ArMBH, 2 )/o'^ 
is the binary’s influence radius and that v a, we 
ca n write the “delay” i n gas-rich nuclear environments 
as ([Granato et al.ll200ll 


^delay, 


gas 


TZffd 


yn 


(29) 


where tdyn = G(Mmbh,i + AfMBH, 2 )/o’^ is the dynamical 
time at the influence radius. In this paper we set TZc = 
10®. This prescription does indeed yield delays < 10^ — 
10® yr in gas-rich environments, as expected. 

In gas-poor enviroments (i.e. Afres < 

ATmbh.i + AfMBH,2), three-body interactions with 
stars do minate and bring the MBHs together on a 
timescale ([Begelman et al.lll980l l 


^delay, stars ^ 5 Gyr 

q 

{l + qY 


4.5 X 10-2 pc 


AImbh,! + All 

hF 


MBH,2 


-1 


, (30) 


where q = AlMBH, 2 /AfMBH,i < 1, and Ogr is the radius 
at which gravitational-wave emission becomes dominant 
over three-body interactions with stars at driving the 
binary’s evolution, i.e. 


Ogr ~ 1.65x10 ^pc X 


(1 + 9)^ 


A/] 


MBH.l 


■AT, 


MBH,2 


(- 


200km/sy GO^AIqPC-® 


lO®Af0 
-il 1/5 


■ (31) 


Here, the density p* of the stellar background in which 
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the binary moves is given by the maximum of the 
stellar-bulge density and the density of the NSC. For 
the latter, we use the average NSC density within its 
half-mass radius, i.e. pnsg ~ [AInsc ~ 2(Mmbh,i + 
AfMBH, 2 )]/[( 8 / 3 ) 7 rr-^], where we decrease the mass of the 
NSC by twice the mass of the binary to approximately 
account for the mass deficit it causes on the NSC during 
its inspiral [cf. equation (l34l) below]. As for the stellar- 
bulge density, our model assumes a Hernquist profile, 
as me ntioned earlier. Com parisons to N-body simula¬ 
tions (jSesana fc Khanll^OlbH show that the appropriate 
radius where this density profile needs to be evaluated in 
order to estimate Ugr correctly is the binary’s mass influ¬ 
ence radius, i.e. the radius at which the enclosed bulge 
mass (in stars) equals twice the binary mass. 

In order to model trip le MBH systems when th ey form, 
we utilize the results of lHoffm an fc Lo^ (l2007f). Bv us¬ 
ing numerical simulations. iHoffman fc Loebl()20(T^ found 
that MBH triples typically trigger the merger of the two 
more massive components on timescales of 10 ® yrs (on 
average) for binary masses of ~ 6 x 1 O®M 0 (and roughly 
comparable masses for the two components), while the 
lightest MBH may be ejected from the galaxy or be left 
wandering far from the galaxy’s center, or less likely come 
back and coalesce with the remnant of the inner binary’s 
merger. Moreover, the time-scales for the merger of the 
inner binary present a log-normal scatter of about 1.4 dex 
around the median value of 10® yrs, as sho wn in Figure 
8 (upper panel) in IHoffman &: LoebI (|2007f) . To extend 
these results to arbitrary masses, we rescale these time- 
scales with the system’s dynamical time at the binary’s 
hardening radius, and we thus obtain an average “delay” 
due to MBH triple interactions given by 


^delay, triple ~ 10 yr 


/ Mmbh,i + AfMBH ,2 ^ 
6 X IO8M0 


(i + g)3' 

(32) 

Whenever a triple forms, we then extract the delay with 
which the two more massive objects merge from a log¬ 
normal distribution centered on this timescale, and with 
r.m.s. of 1.4 dex. For simplicity, we also assume that the 
lightest MBH is ejected from the galaxy and therefore 
lost to the subsequent galaxy evolut ion. 

Finally, IHoffman &: LoelJ (120071) noted that triple 
MBH interactions are less effective at driving MBHs to 
merger in low-mass systems, where the third incoming 
MBH is ejected from the galaxy before it can shrink 
the inner binary to separations at which gravitational- 
wave emission is important. Indeed, the escape veloc¬ 
ity from a galaxy scales as the velocity dispersion, i.e. 


Vesc ^ O' ~ where we defined the mass of the inner 

binary, M = Mmbh.i + A7mbh,2 , and we used the Faber- 
Jackson relation M The separation Ogw at which 

a MBH binary is driven to coalescence by gravitational- 
wave emission within a Hubble time tn can be obtained 
from Eq. (1501) (with the replacements Ogr —>■ Ogw and 


^deiay,stars t Ih): and scales as Ugw ^ ^ Qit (with 

gfj, = MMBH,iAfMBH,2/(AfMBH,i+A7MBH,2)^ the Symmet¬ 
ric mass ratio). In order for the third MBH to be able to 
shrink the inner binary to the separation Ogw, the veloc¬ 
ity V 3 of the third MBH should be lower than Vesc when 
the binary’s separation is Ogw. The energy of the third 


MBH can be estimated simply by energy equipartition as 
E 3 ^ Mmbh. 31'3 GAfMBH,iA^MBH, 2 /agw, from which 
one gets V 3 ~ , with 93 = 

This in turn gives vsjvesc . The nor- 

m alization of this ra t io can be estimated using the results 
of IHoffman fc Lo^ (|2007[) . who find Vesc ~ 1400 km/s 
and V 3 ^ 750 km/s for M 6 x 1 O®M 0 , gj, ^ 1/4 and 
53 ^ 0.25 — 0.3 (cf. their Figs. 5 and 11), thus giving 

3/8 

— ~ 0.5 X - - ^ -- . (33) 

^^esc g 3 /^[M /(6 X 10®)]i/8 

Therefore, whenever this equation gives V 3 > Vesc, we 
assume that the third MBH is ejected without driv¬ 
ing the merger of the inner binary. This implies that 
triple MBH interactions become ineffective at driving 
the merger of binary systems with total mass M < 
2 X 10^ X {ql/ql) Mq. 


Mass deficit caused by black-hole binaries— The 

presence of a MBH binary has profound effects on the 
NSC, causing a “mass-deficit” in the central regions of 
the galaxy. Indeed, during the binary’s inspiral, three- 
body interactions with the stars in the nuclear region re¬ 
sult in an exchange of energy between the binary (which 
becomes more and more bound) and the stars, to which 
large velocities are imparted that are capable of eject¬ 
ing them from the galactic nucleus. Also, when the 
MBH binary finally merges, the resulting MBH rem¬ 
nant acquires a kick velocity up to a few thousands 
km/s due to the anisotropic emission of gravitational 
waves (iCampanelli et al.ll2007ll . and this also contributes 
to remove stars from the galactic core. A similar mass 
deficit may be caused by the ejection of the lightest MBH 
in a triple system. 

We estimate the mass-ejection rate from the NSC due 
to the MBH binary’s inspiral as 

AfJ-' - 

^delay, stars 

(34) 

where q = Mmbh, 2 /A 7 mbh,i < 1 is the binary’s mass ra¬ 
tio and Oh is the the semi-major axis at which the binary 
becomes “hard” (i.e. tightly bound). 

The first term in the numerator of equation (1341) rep¬ 
resents the mass scoured by the MBH binary before it 
becomes hard, where we have identif ied the e j ected mass 
with the mass deficit as defined in iMerrittl (| 200 (tI) ; the 
second term represen ts instead the mass ejected from Oh 
to Ogr (|Merrittll2013f): the explicit expressions for Uh is 
given for instance in lMerritfl (|2013ll . i.e. 


flh- 0.27(1-fg) 


-1 f AfMBH,2\ f 


V IO^A/q J V 200 km/s 


pc , 


(35) 

while for Ogr we utilize equation (1311) . Also, in the de¬ 
nominator of equation (1341) . we note the presence of the 
timescale tdeiay stars, computed via equations (|3ni) and 
m- That timescale accounts for the fact that both 
terms in the numerator where computed for MBH bi¬ 
naries in gas-poor environments, hence the mass deficit 
at the numerator should be “spread” over the timescale 
characterizing stellar interactions. Of course, if the MBH 
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binary lives in a gas-poor enviroment, in our model the 
inspiral lasts exactly tdeiay, stars, and the final mass deficit 
caused by the binary is ^ (Mmbhu + -MMBH,2)[0.7g°'^ -|- 
0.51n(ah/agr)]. If the binary evolution is instead driven 
by gas interactions or by the formation of a MBH triple, 
that mass deficit is suppressed by a factor respectively 

^delay, gas/^delay, stars C>r tdelay, triple/^delay, stars 0 


In addition to the mass deficit caused during the in¬ 
spiral, when the MBH binary finally merges the result¬ 
ing MBH remnant acquires a kick, which, as mentioned 
above, can further remov e mass from the NSC. We es ti- 
mate this mass deficit as (|Gualandris fc Merrittll200^ 


^kick ^ 5 (MmbH,1 + -MmBH, 2) (I4ick/I4sc)^''^^ , (36) 


where I4ick is the recoil velocity of the MBH remnant and 
t4sc is the escape velocity from the galactic center. The 
latter can be easily calculated within our model, from the 
mass and densit y profiles of the bulge and NSC. As for 
I4ick, we f ollow Barausse| ( 201211 and use the analytical 
formula of ivan Meter et ;dT ( 20101) , which fits the results 
of numerical-relativity simulations. 

Similarly, in the case of triple MBH systems, if the 
lightest MBH gets ejected from the system before it can 
cause the inner binary to merge [i.e. , in our model, if the 
ratio vz/vesc given by equation (|M)l is larger than one], 
we assume that the ejected MBH causes a mass deficit 
^kick, triple ^ 5 Mmbh.3, with Mmbh .3 the ejected MBH’s 
mass. We stress that for simplicity we neglect the mass 
deficit caused by the recoil of the lightest MBH in the 
cases in which the triple interactions trigger the merger 
of the inner binary (i.e., in our model, when v^/vesc < 
1). Summarizing, the total mass deficit due to MBH 
coalescences is Mgj = -|- -|- j\^kick,tripie^ 


In addition to the effect of MBH binaries, we also ac¬ 
count for the possible tidal disruption of NSCs by MBHs 
during galaxy mergers. Indeed, if a merger takes place 
between a galaxy (“1”) that contains a NSC (but not 
a MBH) and another galaxy (“2”) hosting a MBH (and 
possibly a NSC), the NSC of “galaxy 1” will be dragged 
by dynamical friction toward the MBH (and NSC if 
present) of galaxy 2, and will therefore be tidally trun¬ 
cated/disrupted, as in the case of star clusters falling 
toward the nucleus of an isolated galaxy. To model this 
effect, we assume that the NSC resulting from such a 
galaxy merger has mass 


A^nsc = AfNSC, 1 X -I- Mnsc, 2 • (37) 

The fraction 0 < < 1 accounts for the tidal trunca¬ 

tion/disruption effects, and is calculated via equations 
m, m, and (f23|l with ax taken to be the velocity dis¬ 
persion of “NSC 1” (which one can compute in terms 
of its mass and size). Note that calculating in the 
same way as for star clusters is justified, at least to first 
approximation. Indeed, as mentioned earlier, the star 
clusters that contribute the most to equation ([2^ have 


^ Note that we have also tried setting the delays between galaxy 
and MBH mergers to very small values tdelay ~ 10^ yr irrespec¬ 
tive of the mechanism driving the binary’s evolution prior to the 
gravitational-wave dominated phase. In this case, the mass deficit 
is always (Mmbh,i + -^mbh,2)[0.7q'^-^ + 0.5 In (ah/ugr)]. This 
test confirms that our results for the NSC evolution are reason¬ 
ably robust, and that the overall conclusions of this paper do not 
depend on our particular model for the delays. 


average initial mass {mgc,in) ~ 2.5 x 10 ®Mq, which is of 
the same order of magnitude as a typical NSC mass. As 
in the case of star clusters, if vt < tk = (with 

c = 30 a typical NSC concentration parameter), we as¬ 
sume that “NSC 1” is fully disrupted, and set = 0. 
On the other hand, if both “galaxy 1” and “galaxy 2” 
have MBHs, we assume that the gravitational field of 
each massive BH protects its NSC from tidal truncation 
and disruption, and we set !F = 1. Similarly, if neither 
galaxy contains a MBH, we set J- = 1. 


3. DATA SAMPLE 

We compare our models to observational data in order 
to gain insights on the reliability of such models. We 
then use our models to put constraints on the origin and 
evolution of NSCs and to understand how their proper¬ 
ties are linked to the evolution of MBHs and their host 
galaxies. 

Our sample of NSC objects was ob t ained by combin¬ 
ing da t a fromiScott fc Grahad (I2013D , lErwin fc Gadottil 
(I2012D . iNeumaver fc Walched (l20fT and by estimating 
NSC obje ct masses of the gal axies in the Fornax cluster 
c atalog ofiTurner et al.l (120121) . 

lErwin fc CadottiT i 20121 ) give mass estimates for a total 
of 18 galaxies that cover Hubble types from SO to Sm with 
dynamically deter mined NSC masses, and 15 Sa and later 
type galaxies from iRossa et al.l (|2006D with NSC mass es¬ 
timates based on high-resolution spectroscopy. Velocity 
dispersi ons for 9 of these galaxies wer e obtained from Ta¬ 
ble 2 of INeumaver fc Walcheil (l2012t) . the rest from the 
Hyper Leda database when available. iScott fc GrahamI 
(I2013D construct ed their sample of NSC s using photo¬ 
metric data from iFerrarese et alJ (l2006aD . iBalcells et al.l 
(|2007D . iCraham fc Snitle d (120001) a nd iCrahamI (|2012bt) . 
Table 1 of IScott fc CrahairJ ( 21)1^ reports 76 galaxies 
containing a NSC with a well determined mass, and total 
mass and velocity dispersion estimates for the majority of 
these. T he uncertainties on the nu c lear object ma s ses are 
give n bv IFerrarese et al.l (l2006aD . IBalcells et al.l l)2007D 
and ICraham fc Soitleil ( 20091 ) as 45%. 33%, and a factor 
of two, respectively. 

We additionally obtained estimates for the stellar 
masses of the full galaxy, Mgx, and the NSC, using the 
photometric data for the 43 early-type galax ies in the 
Fornax cluster catalog of iTurner et al.l (|2012D . To de¬ 


termine, Mgx, we multiplied the total galaxy luminos¬ 
ity in the K band given in the HyperLe da database by 
a ma ss to light ratio of M/Lk = 0.8 (|Bell fc de .lond 
I2001D . based on the typical colors of the bulge pop¬ 
ulation. Galaxy total B-magnitudes, Mb, and yeloc- 
ity dispersions, cr, were also taken from the HyperLeda 
database. The NSC masses were obtained by multiplying 
the published luminosities by the appropriate mass-to- 
light ratio, obtained by using the empirical c orrelations 
betwe en color and mass-to-light-ratio giyen in iBell et ^ 
(I2003D . Typical errors on Mnsc are ~ 35%. 

After eliminating duplicate galaxies contained in more 
than one of the samples - for each galaxy we included in 
the collective sample the mass estimate with the small¬ 
est uncertainty - we collected a total of 89 galaxies with 
reliable estimates for both Mnsc and cr, 127 objects with 
determined Mnsc and Mgx, and 208 galaxies with mea¬ 
sured NSC mass and host galaxy total B-magnitude. 

Although we were not able to obtain all galaxy prop- 











































































12 


Antonini, Barausse and Silk 


cn 

O 


CO 

O 


u 

c/) 


O 


CD 

O 


in 

O 





Fig. 3.— Mnsc against galaxy velocity dispersion (left panel) and galaxy total mass (right panel). Filled circles represent early-type 
galaxies (SO or earlier) while o pen circles are late-type syst ems. Short (long) arrow symbols are NSCs in early (late) type galaxies with only 
an upper limit on their mass JNeumaver &: Walcherir2012l') . Since it is likely that these systems do not contain a NSC we decided to not 
include these upper limits in our regression analysis. Black lines are the computed Loess curves and thin black lines give the corresponding 
95% variance bands. Red lines give the median of the CliN model where only cluster inspirals are considered and corresponding 70% and 
90% confidence-level regions. Blue lines are the results of our semi-analytical model GxeV including cluster inspirals, in-situ star formation 
and processes related to hierarchical galaxy and MBH evolution. Dashed lines indicate the 70% and 90% confidence-level reg ions. Middle 
panel: same as top panel, but with MBH dynamical effects not included in the semi-analytical computations (s ee Se ction 14.1 1 for details) . 
Bottom panel: same as upper panels, but NSCs are allowed to grow only via in-situ star formation (see Section l4.2l for details). 
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Mg. (Msun) 

Fig. 4.— NSC mass fraction as a function of the total stellar mass 
of the host galaxy. Filled circles are early type galaxies, open circles 
are late-type galaxies. Black lines give simple linear fits to the data 
with Mgx < and with Mgx > In the low mass 

regime the NSC to galaxy mass ratio shows a weak correlation with 
Mgx, while in the more massive systems (Mgx > IO^^'^Mq) there 
is a clear trend with Mgx, with the NSC mass fraction decreasing 
rapidly with the galaxy mass. This is consistent with the flattening 
of the scaling relations found at high galaxy masses in Figure [3] 

erties for every object from the literature, our data col¬ 
lectively represent the largest sample to date of NSC 
and host galaxy physical properties. In Figure [3] we plot 
the NSC masses as a function of the velocity dispersion 
of the galaxy spheroid as well as the galaxy mass, for 
early (filled circles) and late (open circles) type systems. 

3.1. Data analysis 

In order to investigate the correlation between NSCs 
and their host galaxy properties we use a locally-weighted 
linear regression m odel (Loess) to ht no n-parametric 
curves to the data (jCleveland 1 1197911198811 . Loess does 
not require one to specify a global function of any form to 
fit a model to the data, as it combines multiple low-order 
polynomial regression models in a fc-nearest-neighbor- 
based metamodel. Unlike a more “standard” simple lin¬ 
ear regression analysis, no assumption is made that the 
data can be htted by a straight line. Loess scatter-plot 
smoothing can therefore be used to reveal complex rela¬ 
tionships that could be overlooked with traditional para¬ 
metric estimation strategies. The obvious trade-off is 
that in general it is not possible to express a Loess model 
with a simple mathematical formula. 

The smoothness of the Loess regression function is de¬ 
termined by the amount of data points used for the lo¬ 
cal fits, a quantity controlled by the span parameter, a. 
Here we select an optimal value of a by using the gen - 
eralized cross validation criterion (GCV: [Graven |[I97^ . 
The basic idea is that the optimal representation of the 
data is obtained by adopting the smoothing pa rameter 
that minimizes the mean-squared error of the fit (jColnb I 

(HiSIlDCMI). 

The solid thick black curves in Figure|3]are Loess inter¬ 
polations obtained with the optimal smoothing parame¬ 
ters a = 0.63 and 0.71 for Mnsc vs cr and Mgx respec¬ 
tively. The 95% variance bands of the Loess curves are 


shown as thin black curves. 

The first noteworthy point is that the Loess curves 
shown in Figure|3]are characterized by a significant bend¬ 
ing at cr > 100 km/s and Mgx > 10^°'®Mq, indicating 
NSC - host galaxy scaling relations that are flat or even 
declining for the most luminous galaxies. Moreover, the 
reconstructed variance bands follow the general trend of 
the Loess curves, indicating that it is unlikely that the 
non-linear features present in the curves are due to ran¬ 
dom fluctuations in the data alone. 

We tested whether or not a linear parametrization of 
the relations provides an adequate description of the data 
by using a T-test t o compare the Loess hts to simpler 
weighted linear fits (|Fox 1119^ . The F tests showed that 
the null hypothesis - i.e., that the bivariate Loess model 
yields no improvement in the fit over the linear regression 
- can be rejected at a high level of significance. For all 
three relations investigated here, the Loess curves pro¬ 
vide a better description to the data than linear mod¬ 
els at > 99% confidence level. This confirms that the 
non-linear features present in the functional dependence 
of Mnsc on the galaxy properties are likely not due to 
noise variability in the data. 

A complementary view is provided in Figure HI where 
we plot the NSC mass fraction as a function of the total 
stellar mass of the host galaxy, i.e., M^sc/Mg^ vs Mgx- 
Evidently, in the low mass regime (Mgx < Mq) 

the NSC to galaxy mass ratio shows a weaker correlation 
with total galaxy mass (with Kendall’s rank correlation 
coefficient tb = —0.28 and associated probability value 
p = 0.06), while in the more massive systems (Mgx > 
Iqio.SMq) there is a clear and steeper trend with Mgx, 
in the sense that the NSC mass fraction decreases rapidly 
wi th the galaxy mas s. This is in contrast with the results 
of IScott fc Graham! (j2013l l. who claimed quite a steep 
decrease of the NSC mass ratio based on fits to the entire 
sample distribution. 

We hnally note that although our analysis demon¬ 
strates that current data are consistent with a signifi¬ 
cant bending of the NSC-host galaxy scaling relations, 
further observations and analysis, combined with theory, 
will be needed in order to refine our results. On the other 
hand, the broadening of the NSC mass distribution for 
the highest luminosity galaxies, as well as the presence 
of low mass nuclei (^ 10 ® Mq) in high mass galaxies 
(-^gx ^ IO^^Mq) certainly poses a serious concern for 
previous claims of the existence of fundamental correla¬ 
tions between NSC and host galaxy properties. More¬ 
over, as we show in the following, the particular form of 
these relations is in agreement with what is expected on 
the basis of our semi-analytical models of NSC formation. 

4. RESULTS 

The red and blue lines in Figure [3] show the re¬ 
sulting NSC mass as a function of a obtained respec¬ 
tively with CUN after 10 Gyr of evolution and with 
GxeV at z = 0. These curves give the median out¬ 
put at a given tr or Mgx as well as its 70% and 90% 
confidence-level regions, i.e. the regions containing re¬ 
spectively 70% and 90% of the NSCs produced by our 
models at a given stellar mass. The CUN model me¬ 
dian can be ht by M^sc ~ 10^Mq((t/50 km/s)^ and 

Mnsc ~ IO^Mq (Mgx/1O^®M0)°'^ at cr < lOOkm/s and 




















14 


Antonini, Barausse and Silk 



I^MBH 

Fig. 5. — Fraction of NSC to MBH mass as a function of MBH 
mass. This plot measures the relative importance of the two types 
of central objects as one proceeds from NSC-dominated to MBH- 
dominated nuclei. The points are observational data and represent 
p^alaxies with measured MBH and NSC mass JGraham fc: Spitleil 
120091 : [Erwin &; Gadottill2013j . The black curve and corresponding 
70% and 90% confidence bands are obtained through our fiducial 
galaxy for mati on model where all the relevant effects described 
in Section 12.21 were included. The blue lines correspond to a 
model in which the scouring effect of MBH binaries forming during 
galaxy mergers was ignored, i.e., we set Mej = 0 in these integra¬ 
tions. Note how the predictions of the two models start to diverge 
at Mmbh ^ 10^Mmq where MBHs and NSCs have comparable 
masses in the model with Mgj = 0; at Mmbh ^ MBH 

mergers become efficient at eroding the surrounding NSC. 


Mgx < respectively. This appears to be in 

good agreement with analytical expectations - 1 Antonlnil 
(|2013f) finds Mnsc ~ IO^Mq (cr/50 km/s)^'^^, adopting 
an idealized isothermal sphere galaxy model. Similar 
scalings and a similar scatter around the median values 
are found with the galaxy evolution model GxeV^ which 
takes into account both cluster inspirals and in-situ star 
formation. The similarity between the results of the two 
models with and without in-situ star formation suggests 
that dynamical friction migration of star clusters is a fun¬ 
damental process for the growth of the nuclei. However, 
as we show below in Section 14.21 dissipative processes 
also play an important role, contributing a significant 
fraction of the total NSC masses in our models. 

Figure [3] shows that at cr > lOOkm/s and Mgx > 
IQio.s^^, scaling relations produced with CUN (i.e., 
a purely dissipationless formation model) appear to flat¬ 
ten toward the most luminous galaxies in agreement with 
the results of the data analysis of Section |3l However, 
a more careful analysis of the plot also shows that such 
a model fails to explain the existence of NSCs of mass 
> IO^AIq which are clearly present in the observational 
data. The reason for this discrepancy is that in CUN 
we have assumed that all galaxies contain a MBH, and 
that these MBHs are in place at the center of all ini¬ 
tial galaxy models before the NSCs grow around them. 
These assumptions are not fully supported by our galaxy 
formation model (see Section 0 and artificially lead to 
an enhanced mass removal from the stellar clusters ac¬ 


creting onto the central NSC, because of the tidal field 
of the MBH. 

The galaxy evolution model CxeV produces scaling re¬ 
lations at z = 0 that are in remarkably good agreement 
with the observed correlations. The scatter around the 
median value clearly increases at the high mass end of the 
NSC distribution. A large population of NSCs are pro¬ 
duced at CT > lOOkm/s that are significantly underweight 
with respect to what the same model would predict by 
simply extrapolating the scaling correlations and scat¬ 
ter from low (Mgx < IO^^ ^Mq) to high galaxy masses. 
In the following section we argue that the formation of 
such “underweight” NSCs and, consequently, the appar¬ 
ent bending of the scaling correlations can be explained 
in terms of the interaction of the NSCs with their host 
galaxy MBHs. 


4.1. The Role of Massive Blaek Holes 

In our semi-analytical models MBHs affect the forma¬ 
tion and evolution of NSCs in two important ways: (i) 
the strong tidal field of the central MBH accelerates the 
mass loss from stellar clusters as they enter its sphere 
of influence and eventually disrupts them as they come 
closer than a distance rdisr, where 


Mmbh 

a’^'disr 


P(0) 


AttG ' 


(38) 


After a cluster is disrupted, most of its stars are dis¬ 
persed around rdisr, which will limit the nuclear den¬ 
sity within this radius. Thus, a preexisting MBH will 
limit the amount of mass that can b e transported t o 
the center by decaying stellar clusters ( AntoniniI20I3H . 
As mentioned earlier [cf. equation p7|) ]. a similar ef¬ 
fect applies to a NSC falling toward a central MBH 
after a galaxy merger. (ii) during galaxy mergers, 
MBH bin aries form and harden by ejec ting surround¬ 
ing stars (jMilosavlievic fc MerrittI l200Ifl . By dynami¬ 
cally heating the NSC, a MBH binary will lower the NSC 
central stellar density, or fully destroy the NSC when 
the bi nary’s mass is significan tly larger than the NSC 
mass (iBekki &: Graharnl I20I0I1 . The kick imparted to 
the merger remnant by asymmetric emission of gravita¬ 
tional wave radiati on will remove additional m ass from 
the galaxy nucleus (|Gualandris fc Merrittll2008ll . 

To demonstrate the role of MBHs in the evolution of 
the nuclei, we performed additional simulations where 
the two dynamical effects due to MBHs mentioned above 
were not included. (Hereafter we refer to these models 
as models “without MBH heating”.) The blue and red 
curves in the middle panels of Figure [3] show the scaling 
correlations resulting from such models. When compared 
to our fiducial models (upper panels), these new integra¬ 
tions produce steeper and significantly tighter Mnsc ~ o’ 
and Mnsc ~ Mg^ relations. This indicates that the break 
in the NSC scaling relations, for which we found evidence 
in the observational data, can be attributed to the inter¬ 
action of NSCs with their companion MBHs. 

Simple arguments can indeed be used to understand 
why the break of the NSC scaling relations occurs at 
a m lOOkm/s. From equation (I38L by requiring rdisr 
xnsc ~ 10 pc, and adopting tk = Ipc and ax = 
20 km/s, we find that MBHs more massive than roughly 
« 10®Mq will suppress the further growth of a NSC via 
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Fig. 6.— The local fraction of nucleated galaxies of our 
models compared to observational constr aints on the nuc leated 
fract ion of early type g alaxies in Vi rgo ([C ote et a . l20Qflh For¬ 
nax JTurner et alJuPlj) . and Coma Jden Bro^^^jh 120^) . The 
lower panels are for late type galaxies. The left panels show the 
results of our fiducial models with all the relevant effects included; 
the right panels correspond to models with no dynamical heating 
due to MBHs (see text for details). A comparison between these 
latter models and our fiducial models indicate that the lack of NSCs 
in galaxies more massive than 10^^ Mq is due to MBHs that 
fully destroy the surrounding clusters during galaxy mergers, and 
also quench their growth by disrupting migrating clusters due to 
their strong central tidal field. 

accretion of stellar clusters. Thus, as the MBH grows, 
the contribution of star cluster inspirals to NSC growth 
decr eases and eventually stops for Mmrh Us¬ 

ing (jFerrarese fc Fordll2005l : Ihbrrarese et al.ll 200 ^ : 


.^MBH ~ 8.6 X 10® 


lOOkm/s 


4.41 




©j 


(39) 


a MBH mass larger than ~ 10®Mq corresponds to stellar 
spheroids with velocity dispersion a > 150 km/s. This 
is consistent with the value of cr at which the NSC-host 
galaxy scaling correlations appear to flatten in the data 
sample. 

Another argument leads us to identify Mmbh ~ 
1 O®M 0 as the critical value of MBH mass above which 
NSCs are significantly affected. The impact of a MBH 
binary on a central cluster depends on the Mnsc/-^mbh 
mass ratio. Given that during a merger a MBH binary 
will eject a mass comparable to its total mass, a require¬ 
ment for a MBH binary to affect significantly the nuclear 
structure of a galaxy is that its mass is of the order of 
or larger than the mass of the surrounding cluster, i.e., 
Mmbh/Mjsisc ^ 1 - 

Using a standard weighted linear regression on our data 
gives the relation tr Ri 50 km/s(MNSc/lO^M 0 )^-®; using 
equation (1391) we find the new relation: 

.^MBH / AInsc 

6 X 1O7M0 ” V6 X IO^Mq/ ■ ^ ’ 



and containing a MBH plus a NSC (dashed lines). The right panel 
shows the fraction of nucleated galaxies, /nd? that contain a central 
AGN, identified here as a MBH with bolometric luminosity L > 
lOio Lq. 

Thus, the mass of the MBH grows faster then the mass 
of the NSC, with the transition from NSC dominated to 
MBH-dominated galaxies occurring at Mmbh ~ 10®M 0 
or (T Ri 150 km/s. This is roughly the value of a at which 
NSCs appear to be significantly affected by their host 
galaxy MBHs in our models. 

The total mass ejected from the center depends also 
on the number of stages in the merger hierarchy that 
have occurred since the MBHs first formed - i.e., after 
Nm mergers, the mass deficit is ~ 0.7 wit h 
Mmbh the hnal (current) black hole mass (lMerrittll2006ll . 
Some low-mass galaxies today might have never experi¬ 
enced a major merger. By contrast, more massive galax¬ 
ies form via mergers of primordial lower mass halos, and 
underwent an above average number of mergers. In ad¬ 
dition, since the occupation fraction of MBHs increases 
with galaxy mass, the scouring effect of MBH binaries is 
enhanced in the higher mass spheroids. 

From FigureO we also see that the GxeV model gener¬ 
ates two distinct populations of NSCs at cr > 100 km/s, 
which can be easily identified by looking at the model 
confidence bands. Nuclei with Mnsc ^ 1 O’^M 0 have 
masses that do not strongly correlate with their host 
galaxy mass and velocity dispersion -- the mass of these 
NSCs lies significantly below the mass that we would ob¬ 
tain by extrapolating the scaling correlations from low 
to high galaxy masses. This population of underweight 
nuclei is a result of the disruptive MBH binary merg¬ 
ers that have partially eroded the surrounding clusters. 
NSCs above the model median seem instead to follow the 
general trend outlined in the low mass galaxy region of 
the plot, where Mnsc ~ cr^- The host galaxies of such 
nuclei did not experience a major merger event since the 
epoch of the last major gas accretion event, so that their 
NSCs remained essentially unaffected by MBH binaries 
up to the present epoch. This picture is in agreement 
with expectations based on hierarchical models of galaxy 
formation - for example, iFakhouri et al.l (120101) find that 
for a dark halo mass of 10^^ Mq only 31, 53 and 69 per 
cent of these halos have experienced a major merger since 
z =1, 2 and 3, respectively. 

We can now ask which of the two processes, tidal dis¬ 
ruption of star clusters or mass ejection from MBH bi¬ 
naries, is responsible for the broadening of the scaling 
correlations found in the data. 
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We find that both CUN models with Mej = 0 and 
ones without the tidal truncation/disruption of clusters 
by the central MBH produce scaling correlations that 
appear to flatten at cr > 100 km/s. The NSC scaling 
correlations produced by GxeV models in which we set 
Mej = 0 (but in which we included the tidal disruption 
of migrating clusters) are instead at odds with obser¬ 
vations, as they show no broadening of the NSC host 
galaxy property correlations for high velocity dispersion 
galaxies; only GxeV models that do take into account 
the scouring effect of MBH binaries were found to be in 
good agreement with the observed relations. We con¬ 
clude that in GxeV the crucial ingredient to reproduce 
the MNSc-host galaxy property scaling relations is the 
scouring effect of MBH binaries, while the tidal strip¬ 
ping of migrating clusters by central MBHs is only a 
secondary effect in these models. We remark that the 
GxeV model has a clear advantage over the CUN model, 
because the former follows the hierarchical evolution of 
NSCs and MBHs along merger-trees, while the latter as¬ 
sumes a monolithic evolution. Although the results of the 
two models are generally in good agreement with one an¬ 
other, the differences outlined above are mainly a result 
of the implicit assumption made in the CUN model that 
the NSCs always grow around pre-existing MBHs. This 
leads to an artihcially enhanced mass removal from the 
stellar clusters accreting onto the central NSC, because 
of the MBH tidal field. 

To illustrate more clearly the role of MBH binaries in 
shaping their companion NSC properties, we compare in 
Figure [S] the observationally constrained NSC to MBH 
mass ratio in galaxies containing both types of central ob¬ 
jects to the same ratio obtained (i) in our fiducial GxeV 
model containing all the relevant effects described in Sec¬ 
tion [521 and (ii) a model in which we set Mej = 0, i.e., 
in which the scouring effect of MBH binaries was not 
included in the calculation. Clearly, the scouring effect 
due to MBH binaries described by equations (IM)) and 
(l36l) is the crucial ingredient to reproduce the observed 
correlation. Also, note that the correlations produced 
by the two models start to diverge from each other at 
Mm rh ~ lO^Mfn, in agre ement with our predictions (see 
also I Antonini et al.l 1201511 . 

4.1.1. NSC and MBH occupation fraction 

Observationally, the frequency of nucleation in early 
type galaxies is found to increase sharply from zero for 
spheroids brighter than Mb = —19. 5 to > 90% for 
galaxies fainter th an this magnitude (jCote et al.l 120061 : 
[Turner e t al.ll20 12ll. In a sam ple of 332 late type galax- 
ies. iGeorgiev fc BOkerl (|2014ll found that > 80% of these 
galaxies harbor a well defined NSC. Hence NSCs are 
found in most galaxies of all Hubble type, but tend to 
disappear in the brightest spheroids. 

Figure|n]compares the frequency of nucleation obtained 
in our models to that of late and early-type galaxies, 
as derived from observations. We define here early- 
type galaxies as systems with bulge-to-total mass ratio 
Mbuige/Mgx >0.7, and late-type galaxies as systems with 
Mbuige/Mgx ^ 0.7. In agreement with observations, our 
models predict that almost all early-type galaxies that 
are less massive than Ri 10^^ Mq contain a NSC, and 
that the frequency of nucleation is nearly zero for galax¬ 
ies with mass larger than this value. We note that the 


exact occupation number obtained through GxeV should 
be considered somewhat approximate - the identification 
of a galaxy as nucleated in this model is uncertain given 
that we cannot directly reconstruct the density profile 
of the NSC and compare it to the density profile of the 
galaxy background. Nevertheless, the occupation frac¬ 
tion is found to decrease significantly for early type galax¬ 
ies with mass larger than 10 ^^ Mq, in good agreement 
with observational constraints. The bottom panel of Fig¬ 
ure [5] shows the frequency of nucleation in late type sys¬ 
tems. Approximately ^ 90% of these galaxies were found 
to contain a NSC, also in fairly good agreement with ob¬ 
servations. 

The observed drop in the nucleation frequency at 
high galaxy masses is often attributed to the disrup¬ 
tive effects of MBH binaries forming during galaxy merg¬ 
ers (IBekki fc: Graham~ll2010ll and to the tidal disruption 
of migrating cl usters at large galactocentric distances by 
central MBHs (|Antoninill2013il . In the right panels of Fig- 
ure|6]we show the nucleated fraction in models where the 
dynamical heating due to MBHs on migrating clusters 
and during mergers was not included. In these models 
the percentage of nucleation remains close to 100% re¬ 
gardless of galaxy luminosity, demonstrating that MBHs 
are responsible for the absence of NSGs in the most lu¬ 
minous galaxies, as also argued in previous work. 

NSGs and MBHs are known to coexist in some 
galaxi es across a wide range of masses and Hubble 
types (ISeth et al.ll2008all . However, not all NSGs contain 
a MBH; an example is the M33 NSG that has no central 
MBH down to highly const raining limits (iMerritt et al.l 
I200II: iGebhardt et al.! 1200111 . Studies indicate that at 
least some NSGs can host a central MBH, but the overall 
occupation fraction of MBHs in NSCs remains unknown. 
The left panel of Figure |7| displays the fraction of all 
galaxies in the GxeV model that contain a MBH and the 
fraction of galaxies containing both a MBH and a NSC. 
Our models predict that a high fraction of galaxies of 
intermediate and low luminosity contain a central NSC, 
but that the NSCs tend to disappear in massive galax¬ 
ies. Thus, the total fraction of mixed systems with both a 
NSC and a MBH is roughly equal to the MBH occupation 
fraction for galaxies with mass Mgx < 10^^ Mq, and is 
roughly equal to the NSC occupation fraction for galax¬ 
ies more massive than 10^^ Mq. In the right panel of 
Figure |7| we show the fraction of nucleated galaxies that 
contain a MBH with bolometric luminosity L > 10^°Lq. 
We find that the fraction of galaxies with a NSC that also 
have an active MBH increases from ~ 5% to 30% going 
from Mgx ^10® Mq to ^ 10^^ Mq. These results are 
in fa irly good agreemen t with the observational results 
from ISeth et al.l (I2008all . who found that about 10% of 
NSCs in their spectroscopic sample also host an AGN. 

4.2. In situ vs migration: the relative contribution to 
NSC growth 

The lower panels of Figure [3] show the scaling rela¬ 
tions generated by a version of the GxeV model with 
no contribution from cluster inspirals. In these models, 
the only mechanism responsible for the formation and 
growth of the nuclei is local fragmentation of gas that 
leads to star formation and in-situ build up of a central 
compact cluster. 

The lower panels of Figure [3] demonstrate that even a 
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Fig. 8. — Ratio of the median mass obtained in the only in-situ 
star formation model to the median mass of our fiducial model 
which also includes cluster inspiral processes ((Min—situ)/(-^NSc))- 
Error bars represent Icr uncertainties. A large fraction of the mass 
of our model NSCs is built up locally from episodes of star forma¬ 
tion triggered by infalling gas. 



purely dissipative model provides a quite good descrip¬ 
tion of the observed NSC-host galaxy scaling relations. 
The median of the mass distribution as a function of cr 
can be fit by Mnsc ~ 6.5 x 10 ®(ct/ 100)^'^M©. Hence, 
in this model we obtain scaling relations that are consis¬ 
tent with the observation that NSCs follow scaling rela¬ 
tions that are shallower than the corresponding ones for 
MBHs. 

We can derive analytical scaling relations for the in- 
situ formation model by considering the MBH and reser¬ 
voir characteristic timescales of evolution. Let us con¬ 
sider a situation in which the MBH is hosted by a gas- 
rich nucleus [resulting from a recent starburst, cf. equa¬ 
tion (IMl) ]. i.e. Mmbh ^ Afres- The MBH will then 
grow approximately at the Eddington rate, i.e. Mmbh ~ 
MMBn/tsaip, with tsaip the Salpeter timescale 

= £5 = (ot) 

where e is the radiative efficiency and k = 0.398cm^ g“^ 
the electron scattering opacity. Gas accumulates at the 
center of the galaxy mainly during starbursts [cf. equa¬ 
tion (l24l) ] . Bursts of star formation will take place on the 
dynamical timescale of the galactic bulge, i.e. 

(TSoSvi) 

with R the galactic bulge scale radius. In the last expres¬ 
sion we have used the fact that both R and cr are related 
to the tot al galaxy luminosity through the empiric al cor¬ 
relations (jBernardi et al.l[200^ iDekel fc Coxll2006f) : 

(«) 

and 

/ L \ 



Fig. 9. — Scaling correlations of NSCs and MBHs in our galaxy 
formation model at different redshifts. At high redshift the NSCs 
are the dominant central component of galaxies; between 2 = 3 
and 2 = 1 the MBHs grow faster and by 2 = 2 they are the 
dominant central component in galaxies with cf > 100 km/s. After 
this point, MBH mergers become efficient at carving out the pre¬ 
existing nuclei. This effect induces the bending/broadening of the 
relation toward high galaxy masses at low redshifts, which 
can be clearly seen in the bottom panels of the figure. In the 
bottom right blue points repr esent NSCs and black points represent 
MBHs lITremaine et al.|[2002l) 

We can thus write, from equation (|24)) . Mres ~ 
HresM®i[jgg/tdyn during the starburst event, and there¬ 
fore 

AMres ^ Mres ^ ^salp ^res 

AMmbh Mmbh Wn Mmbh 

where AM^es is the total mass of cold gas that falls to 
the nucleus during the star formation event, and AMmbh 
and AM^^jgg denote the changes in the MBH and stellar 
bulge masses during the starburst. (Note that eventu¬ 
ally all the gas of the bulge is transformed into stars, if 
feedback is neglected, i.e. M®“gg Ri AM^^j^^.) 

Now, if we assume that the MBH approximately 
satisfies the Magorr ian relation (iMaeorrian et al.l 119981 
iHaring fc Rixl |2004) Mmbh ~ [with Am Ri 
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Fig. 10. — Mass distribution of MBHs and NSCs at various red- 
shifts. At high redshift, 2 : < 4, the NSCs are typically more mas¬ 
sive than MBHs. At 2 : = 0 NSC masses are typically of the or¬ 
der with only a few rare clusters having masses above 

10® Mq. Note the transition from NSC dominated to MBH domi¬ 
nated galaxies for Mmbh ^ 10®Mq occurring at 2 : < 2. 


1.4 X 10-3 (iHaring fc Riyll2nnl ]. we can write 


^-^res ^ tsalp 

AMmBH idyn 




salp 


Ar. 


^dyn A 


M 


5 r 

"/crit OC (7 


-2.06 


(46) 


where have assumed AM,^^jgg/M^^jgg Ri /crit, because 
starbursts are typically triggered by major galaxy merg¬ 
ers, i.e. ones with mass ratio between the baryonic com¬ 
ponents of the two galaxies larger than /crit ~ 0.25. 0 
Eventually all the gas of the reservoir will either form 
stars or accrete onto the MBH, i.e. AMres = AMmbh + 
AMnsc- However, since we assumed Mmbh “C Mres, 
most of the nuclear gas will end up in the NSC rather 
than in the MBH, hence AMres ~ AMnsc, and we can 
write, for Mmbh < Mnsc, 


We can then integrate this equation to give 


Mnsc 


1.9 X ^4:31Mmbh 


t, 


dyn 

7.8 X IO^Mq 


lOOkm/s 


2.35 


(48) 


where we have used equation dMl) and assumed a ra¬ 
diative accretion efficiency e ~ 0.1. Although very sim¬ 
plified, this model approximately reproduces the scaling 
relations generated by our dissipative GxeV model (for 
Mmbh ^ Mnsc) in slope, and also (within a factor of a 
few) in normalization. 

The mass growth of NSCs is regulated by two pro¬ 
cesses: inspiral of star clusters and in-situ star forma¬ 
tion. In order to quantify the amount of NSC mass that 
is contributed by local star formation in our models, we 
computed the ratio of the median NSC mass obtained by 
including both cluster inspirals and in-situ star formation 
to that obtained with only the contribution of in-situ star 
formation. We plot this quantity as a function of total 
galaxy mass in Figure El 

In our models, about half of the mass of NSCs is con¬ 
tributed by local star formation for galaxies less massive 
than ~ 3 X IO^^Mq, while for galaxies more massive than 
this value, the contribution from dissipative processes be¬ 
comes progressively more important. Some provisional 
evidence for a preferential dissipative mode in galaxies of 
progre ssive ly larger mass e s was presented bv iCote et al.l 
i)2006f l and iTurner et al.l 1)2012( 1. These authors showed 
that the brightest nuclei which reside in the most lumi¬ 
nous hosts have colors that are significantly redder than 
expected on the basis of a star cluster merger scenario, 
suggesting that an additional process - e.g., dissipational 
infall of metal-rich gas during mergers - likely begins to 
dominate the formation of nuclei in galaxies of higher 
masses. 

We finally note that the relative contribution of the 
two formation channels depends on the adopted value of 
star cluster formation efficiency /gc which remains a quite 
uncertain parameter of our models. However, given that 
plausible values of /gc were found to only impact the 
normalization of our results for the empirical relations 
and to not affect their slope, assuming /gc ~ const we 
can simply express the relative contribution of the two 
formation channels using the general formulation: 


dMNsc 

dtr 


^salp A,.es 

^dyn Am 


/crit 


dMMBH 

dtr 


(47) 


(Mci) 

(Min —situ) 


1.5 X 



(49) 


where we used AMnsc ~ (dMNsc/dcr)ACT and 
AMmbh ~ (dMMBH /dcr)A(70 


with (Min-situ) the typical NSC mass that originates in- 
situ from episodes of star formation and (Md) the mass 
brought in by migrating clusters. 


Note that the same approximate scaling AMres/AM mbh oc 
cr”2 06 obtained from equation ll45l l by assuming that the 

MBH satisfies the M — a relation Mmbh if observes 

that (^bulge/^)^ (where we have used equation 

(i44t and assumed constant mass-to-light ratio). Indeed, this latter 
relation yields ~ ct^Act, which replaced in equation 

gives AMres/AM mbh oc if one assumes Act ct in major 

galactic mergers. 

^ Note that dMMBn/dcr and dMNSc/<i<^ essentially the rates 
of change of the MBH and NSC masses during a starburst (time 
being parametrized by the host galaxy’s velocity dispersion a). 


5. COSMOLOGICAL EVOLUTION OF SCALING 
CORRELATIONS 

The existence of a fundamental con nection between 
MBHs and NSCs was first suggested bv iFerrarese et all 
(|2006all . These authors showed that NSCs and MBHs 
follow similar scaling relations with their host galaxy 
properties, and thus argued that they are two differ¬ 
ent manifestations of the same astrophysical type of sys- 
tem, which the y called “central massive object”. Later, 
iGrahaml ()2012f) used a large sample of nucleated galaxies 
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to show that the scaling relations of NSCs and MBHs are 
quite different from each other, with the former follow¬ 
ing much shallower correlations with their host galaxy 
properties. Graham’s hndings might suggest that NSCs 
and MBHs formed through different physical processes. 
This could be the case if for example NSCs originated 
elsewhere in the galaxy and then migrated to the cen¬ 
ter th rough dynamical friction processes fe.g.. lAntoninil 
[ 2 (^. _ 

iNenmaver fc Walcheil (|2012f ) presented a first Mmbh 
versus Mnsc diagram and found, in agreement with our 
study, the existence of three different regimes: (a) NSC 
dominated nuclei at tr < 100 km/s, (b) a transition re- 
gion, and (c) MBH-dominat ed nuclei at cr > 150 km/s. 
iNeumaver fc Walched (I2012D argued that this is consis¬ 
tent with a picture in which black holes form inside NSCs 
with a low-mass fraction. They subsequently grow much 
faster than the NSC, destroying it completely when the 
r a.tio grows abo ve ^ 100. 

INeumaver fc Walcherl (|2012fl also argued against MBH 
mergers as responsible for the disruption of the nuclei 
in the highest mass galaxies. These authors pointed 
out that if NSCs are disrupted during mergers, elliptical 
galaxies - thought to be the product of galaxy mergers 
- should rarely host a NSC. The fact that most early- 
type galaxies have a NSC would therefore suggest that 
mergers do not play a pivotal role in leading to NSC 
disruption. We do not agree with this interpretation. 

As discussed in Section 14.11 a MBH binary will 
eject from the galaxy center a mass comparable to its 
own mass, so that only mergers with MBH to NSC 
mass ratio larger than unity can significantly heat the 
NSC and make it susceptible to destruction during the 
merger event. NSCs can therefore survive and indeed 
grow during the morphological transformation from disk- 
dominated to bulge-dominated galaxies if the progenitor 
galaxy NSCs are significantly more massive than their 
central MBHs. Indeed, from the NSC and MBH scaling 
correlations, we hnd that NSCs dominate the nuclei of 
galaxies with velocity dispersion a < 150km/s, so that 
mergers of galaxies at the low end of the a distribution 
will not lead to the complete disruption of the host galaxy 
NSCs, while mergers of massive galaxies, characterized 
by a larger MBH to NSC mass ratio, will lead to their 
full disruption. This picture is consistent with both the 
abundance of NSCs in early type galaxies of intermediate 
luminosity, and with the fact that NSCs tend to disap¬ 
pear in the brightest spheroids hosting the most massive 
MBHs. 

To illustrate the simultaneous evolution of MBHs and 
NSCs in our semi-analytical galaxy formation models, 
we plot in Figure [9] the MBH and NSC scaling correla¬ 
tions at different redshifts and the corresponding mass 
distributions in Figure ITOl Specihcally, in order to calcu¬ 
late the median and confidence regions in our model we 
have have only considered the MBHs residing in bulge- 
dominated galaxies (which we identify with ones having 
bulge to total mass ratio larger than 0.7)0 

We find that at high redshifts, 2 > 3, the nuclei of 
galaxies are dominated by NSCs. Between redshift 4 


® We make this choice to correct for the observational bias that 
selects galaxies with a significant bulge component, so that mea¬ 
surements of a are possible in the first place. 
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Fig. 11.— Same as Figurel^ but for the heavy MBH seed scenario 
of lKoushiappas et al.l ll 200411 . 

and 2 the MBHs grow faster than the NSCs, becoming 
by z Ri 3 the dominant nuclear component in galaxies 
with a > lOOkm/s. After this point, the NSC scaling 
correlations start to flatten at high values of a, as MBH 
binaries forming during mergers are now efficient at erod¬ 
ing the surrounding clusters. In addition, the NSCs can 
no longer grow efficiently in the most massive galaxies, 
since inspiraling stellar clusters are tidally disrupted at 
larger galactocentric distances in galaxies with progres¬ 
sively more massive MBHs. 

Figure [TO] displays the redshift evolution of the mass 
distribution of MBHs and NSCs. In these plots we in¬ 
clude all galaxies in our models regardless of their spe¬ 
cific value of bulge to total mass ratio. At high red¬ 
shift the NSCs are typically more massive than MBHs; 
after z r; 2 the mass distribution of NSCs intersects 
at about 10® Mq the distribution of MBHs. This sets 
the transition between NSC dominated and MBH dom¬ 
inated galaxies seen also in the observational data. We 
also compared our synthetic NSC mass distributions for 
early type galaxies with the observed mass dis tribution of 
NSCs for t he catalogs of early ty pe galaxies of iCote et al.l 
(|2006f) and iTurner et all (|201^ . Anderson-Darling and 
Kolmogorov-Smirnov statistical tests used to compare 
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Fig. 1 2.— Same as Figur e I10| but for the heavy MBH seed sce¬ 
nario of IKoushiappas et al.l II2004I '). 

the observed to the MBH model cumulative distributions 
give p-values of 0.35 and 0.42 respectively, indicating that 
the hypothesis that the simulated and observed distribu¬ 
tions are significantly different from each other can be 
rejected at a high level of confidence. For a compari¬ 
son wit h the observed M BH mass function at z = 0, see 
instead iBarauss^ 1)20121 1 and iSesana et al.l l)2014r i . 

5.1. Dependence on the MBH seed model 

In the previous sections, we have discussed the results 
of a galaxy formation model in which MBHs are evolved 
starting at z ^ 15—20 from initial “light” masses Mgeed ~ 
200 Mq. This could be the case if MBHs fo rmed as 
remnants of Pop HI stars ()Madau fc Reesl[2001ll . 

In order to test the dependence of our results on the 
assumed scenario for the formation of MBHs, we evolved 
additional models in which halos are initially populated 
by MBH seeds of mass Mgeed ~ 10® Mq, which could for 
example form as a result of the collap se of massive proto- 
galac t ic disks at redshifts z 10 — 15 (IKoushiappas et ahl 
2004 iBegelman et al.l l2006riLodato fc NataraianI 120061 
Volonteri et al.ll2008ll . In particular, we consider a model 
in which these “h eavy” seeds form with h igh halo oc¬ 
cupation number (IKoushiappas et aT1l2004l . and one in 
which the hal o occupation number i s moderate (namely 
the model of iVolonteri et al.l (120081) , where we set the 
critical Toomre parameter below which the formation of 
a MBH seed becomes possible to Qc = 2). 

Figured!] shows the redshift evolution of the Mnsc — cr 
and Mmbh — u r elations in the case of the MBH heavy 
seed scenario of IKoushiappas et al.1 (|2004ll . Figure 
displays the corresponding MBH and NSC mass distri¬ 
butions from z = 4 to z = 0. A comparison of these 
plots with those in Figures |9| and [101 demonstrates that 
the local mass distribution of NSCs and their evolution 
from high redshift is not affected by the MBH seed model 
in any important way. (Of course, the MBH mass func¬ 
tion in the heavy-seed models differs from that in the 
light-seed model at the low-mass end, which is domi¬ 
nated by MBHs that have not evolved significantly from 
their seeds. However, at intermediate and high masses. 


the MBH mass functions are very weakly dependent on 
the seed model, since memory of the initial conditions 
has been lost due to accretion and mergers.) Note that 
also in the heavy-seed scenarios, we find that the NSC 
cumulative mass distributions at z = 0 for early type 
galaxie s are consis t ent w ith th e observed m a ss dis tribu- 
tion in iCote et ^ (1200611 and [Turner et al.1 (1201211 at a 
high level of si gnificance (simi l ar con clusions hold for the 
seed model of Volonteri et al.1 (|2008ll l. 

The fact that our results about NSC evolution are not 
sensitive to the particular model chosen for the forma¬ 
tion of MBH seeds is primarily a consequence of the high 
NSC to MBH mass ratio at high redshift. At z > 3, 
in both the heavy and light seed scenarios, the NSCs are 
the dominant central component of galaxies (with the ex¬ 
ception of small NSC masses < 10 ®Mq, at which MBHs 
may dominate if they are seeded with high halo occupa¬ 
tion number at high redshift). It is only after the peak of 
the quasar activity at z ~ 2 that MBHs become massive 
enough to significantly affect the NSCs at intermediate 
and high masses. However, by this time the MBH mass 
distributions in different seed scenarios are very similar 
to each other (again, with the exception of the low-mass 
end), and therefore the subsequent evolution of the NSCs 
is also very similar. The mass growth of the MBH popu¬ 
lation is in fact dominated by the mass accreted during 
the quasar epoch at z Ri 2, thus washing out the imprint 
of the initial conditions, cf. also iBaraussd ()20I2f) . 

Since to date the formation process of MBHs remains 
largely unconstrained, the fact that our results are ro¬ 
bust against different initial seed formation scenarios is 
important if we want to make robust predictions about 
the overall evolution of the NSC population. 

6. DISCUSSION 

6 .1. Formation of the central regions of galaxies 

In this paper, we have presented a study aimed at un¬ 
derstanding how the central regions of galaxies formed, 
and how the evolution of MBHs and NSCs is connected 
to that of their host galaxy. 

In FiguredSl we show examples of NSC and MBH evo¬ 
lution predicted by the GxeV model. These cases should 
be seen as qualitative because of the great variety of pos¬ 
sible NSC and MBH histories that are possible within our 
model, but we have chosen examples that are roughly 
representative of the various possibilities mentioned be¬ 
low. In more detail, we consider the evolution of NSCs 
back in time following their main progenitor, i.e. we start 
from the NSC hosted in the central galaxy at z = 0, and 
at each galaxy merger we follow its more massive NSC 
progenitor. When no NSC is present at z = 0 in the cen¬ 
tral galaxy (namely in the bottom-left panel), we follow 
the main-progenitor history of the MBH. As can be seen, 
NSCs grow by the combined action of (i) stellar-cluster 
infall, which happens continuously throughout their his¬ 
tory, but is particularly enhanced in starburst galaxies 
(cf. the blue line in Figure |T21 which denotes the mass 
in stellar clusters, and which shows a smooth evolution 
with superimposed spikes due to starbursts); and (ii) in- 
situ star formation, which takes places for the most part 
in starburst galaxies formed by major mergers (cf. the 
spikes in the mass of the low-angular momentum reser¬ 
voir available for nuclear star formation, i.e. the green 
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Fig. 13. — Examples of possible main-progenitor evolutions of NSCs and MBHs, as predicted by our model and for different final products 
at 2 = 0. At upper left we indicate the final product of the evolution, a composite NSC H-MBH nucleus (upper panels), a NSC without 
MBH (middle panels), or a MBH dominated nucleus lacking a NSC (bottom panels). We consider the evolution of the NSC hosted in the 
central galaxy at z = 0, and proceed back in time, following the main NSC progenitor at each merger. However, when no NSC is present 
in the central galaxy at 2 = 0 (namely in the bottom-left panel), we follow the MBH main-progenitor history. 


line in Figure [T51) . Note also that MBH mergers have a 
prominent effect on the MBH and NSC history, not only 
because MBH binaries erode or even completely destroy 
the NSC if the MBHs are sufficiently massive, but also 
because the merger can eject the remnant MBH from the 
galactic nucleus, when the kick velocity imparted by the 
anisotropic gravitational-wave emission is larger than the 
escape velocity from the central parts of the galaxy. Note 
instead that the apparent “re-growth” of the MBH after 
an ejection, which can be seen in Figure 1131 is simply 
a consequence of our choice of following the NSC main- 
progenitor history: at certain galaxy mergers, the galaxy 
containing the main NSC (but no MBH) may merge with 
one carrying a MBH. 

Below we discuss in more detail how the evolution of 
the galaxy, and its merger history can give rise to the 
variety of galactic central structures observed today. 

Nucleated galaxies containing a MBH (cf., upper pan¬ 
els of FiaureUM. 

If the MBH mass does not grow above ~ 10^ Mq 
throughout the evolution of the galaxy, the central NSC 
cannot be significantly eroded during galaxy mergers by 
inspiraling MBH binaries, because of the large NSC to 
MBH mass ratio. Moreover, stellar clusters will migrate 
to the center essentially undisturbed by the tidal field of 
the central MBH. 

If during the last merger event the MBH remnant is 
retained, the end product will be a galaxy containing 
both a central NSC and a MBH; An example of such a 


galaxy is the Milky Way which hosts a MBH of mass 
ft! 4 X 10® Mq (jChez et al.l 119981 iCillessen et al.l 1^00911 
and a NSC of nrass ft 3 x 10^ Mq (|Sch5del et al.l[20M 
iFeldmeier et al.l 12014( 1 . A handful of external galaxies 
are also know to contain both a NSC and a MBH, which 
are oft en found to have comparable masses (jSeth et al.l 
l2008bf) . These galaxies lie near the transition region be¬ 
tween MBH dominated and NSC dominated systems. As 
the MBHs of these galaxies grow, MBH mergers will par¬ 
tially destroy their host NSCs in the process (e.g., upper 
right panel of Figure fl^ . 

Nucleated galaxies without MBH (cf., middle panels of 
Fiaure[TS\}. 

After a merger a MBH can receive a “kick” due to gravi¬ 
tational recoil with velocities as large as a few thousands 
km/s, which may eject the MBH from its host galaxy. 
The ejected MBH will carry a cluster of bound stars, 
a hypercompact stellar system that can appear similar 
in size and lu minosity to star clus ters or ultracompact 
dwarf galaxies (iMerritt et al.l[200^ . If the ejected MBH 
has a mass < IO^Mq, the NSC will remain virtually 
unaffected. If from that point on the galaxy evolves pas¬ 
sively without experiencing mergers which might bring a 
new MBH to the center, the evolution’s end product will 
be a nucleus containing a NSC but no MBH. Examples of 
such type of evolution are illustrated in the middle panels 
of Figure [T51 A similar end product can be attained if a 
MBH seed never formed and the galaxy evolved without 
experiencing major mergers during its entire evolution. 
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We note in passing that whether a MBH seed will form 
might depend s on the properties of t he central cluster at 
high redshift (iMiller fc Daviesl2012ll . If the accumulated 
cluster has a high velocity dispersion cr > 40 km/s, this 
might seed the growth of a central MBH because kine¬ 
matic heating from binary stars is insufficient to prevent 
complete core collapse, which might then lead to run¬ 
away mergers and the formation of a central MBH seed 
of ^ 10^ Mq. In galaxies hosting clusters with initial 
velocity dispersion less than 40 km/s, binaries will pre¬ 
vent core collapse and consequently the formation of a 
MBH seed. Galaxies with clusters with initial velocity 
dispersions below this limit might therefore never form a 
MBH. Such evolutionary path might lead to the forma¬ 
tion of galactic nuclei such as those of M33 or NGC205, 
which lack a MBH and host a low velocity d ispersion 
NSG (jMerritt et al.ll200'll : iGebhardt et al.ll2001h . 

MBH dominated galaxies without NSC (ef., lower pan¬ 
els of FiaureYWi). 

After the MBH mass grows above > 10® Mq, any ma¬ 
jor merger will be highly disruptive for a central NSC; 
MBH binaries forming during major mergers in this high 
MBH mass regime will carve out the galactic center of 
stars destroying a pre-exi sting NSC, eventually produc¬ 
ing a central mass-defic it (|Milosavlievic fc Merrittl[2001l : 
iBekki fc Grahain1l20intl . 

The accretion of a nucleated dwarf galaxy by a low- 
density giant galaxy, could brin g a NSC to the center 
of the latter galaxy. However, iMerritt fc Crud 1)20011) 
showed that the secondary galaxy is disrupted during the 
merger by the giant galaxy MBH tidal field, producing 
a remnant with a central density that is only slightly 
higher than that of the gian t galaxy initially. Moreover, 
as shown in I Antonini (120131) . after a NSC is disrupted in 
a dry merger it will have a hard time to regrow, because 
the MBH tidal field is effective at disrupting migrating 
clusters. Thus, once the galactic MBH mass grows above 
> 10® Mq, the central NSC is likely to be disrupted 
and the central core is likely to be preserved during the 
subsequent evolution of the galaxy. 

Low mass spheroids without NSC. 

NSCs tend to disappear in Galaxies fainter than My ^ 
— 12 (jvan den Bergh l|l98fil) . While the purely dissipa¬ 
tionless formation model CUN predicts that the NSC 
occupation fraction decreases substantially in low mass 
systems in agreement with observations, the GxeV model 
results in a high occupation fraction of NSCs in low mass 
galaxies, and it is therefore in tension with observational 
findings. We believe that the reason for this discrep¬ 
ancy is that discreteness effects are not accounted for in 
GxeV as argued in the following. In this latter model, 
we compute the flux of mass accreting onto the nucleus 
by averaging over the initial mass function of the stellar 
cluster population. In GUN, we simulate instead the in¬ 
spiral of each of the formed clusters, drawing their masses 
from the same initial mass function. Because of the ini¬ 
tial small number of stellar clusters in low mass galaxies, 
and because of the steep power law that we adopt for 
their initial mass function, some of these galaxies will 
host no cluster that is massive enough to decay to the 
center before being disrupted in the process. This sug¬ 
gests that the lack of NSCs in spheroids less luminous 
than My ~ —12 is a consequence of the initial low num¬ 


ber of massive clusters in these galaxies. Clearly, this ef¬ 
fect cannot be accounted for in the GxeV model, which 
does not follow the evolution of individual stellar clusters 
but rather the average mass infall to the nucleus. 

6 .2. Galaxy cores and eusp regeneration 

Our analysis focuses on the mass evolution of NSCs 
and ignores effects due to dynamical relaxation that can 
change the density profile of the clusters over time. We 
have shown that the merger history of galaxies after z ks 
2 plays a key role in determining the structure of NSCs 
observed today and their disappearance in the brightest 
galaxies. In our models, we did not account for relaxation 
processes that can rebuild an eroded cusp as stars diffuse 
into the MBH sphere of influence. Even if a NSC is 
fully disrupted during a merger, a stellar cusp of mass of 
order 10% the mass of the black hole and extending out 
to one tenth its sphere of influence wi ll reform over the 
relaxation time (IMerritt fc Szelll I2006D . Although cusp 
regrowth is not expected to have an important effect on 
the mass evolution of the NSCs, it is worth discussing it, 
as the ability of a MBH to shape a galactic nucleus and 
the observational consequences depend on whether cores 
a re regenerated or p r eserve d after mergers. 

iBekki fc Grahami ()2010ll performed A^-body simula¬ 
tions to investigate the evolution of merging star clus¬ 
ters with seed black holes. Using black holes masses of 
1 — 5% their host NSC ones, they investigated how the 
binary coalescence dynami cally heats and de s troys the 
surrounding clusters. While lBekki fc Graham I (|2010ll ar¬ 
gued that this mechanism provides a pathway to explain 
the observed reduction in the nucleus-to-galaxy stellar 
mass ratio as one proceeds from dwarf to giant elliptical 
galaxies, they did not discuss how their results would be 
affected by two-body relaxation after the MBH binary 
merges, and between mergers. 

During galaxy mergers, cores are carved out in the 
stellar distribution, with radii of order the influence ra¬ 
dius of the massive binary. This makes NSC suscep¬ 
tible to disruption. On the other hand, if the supply 
of stars is continuously replenished, the damage to the 
NSC can be much smaller. Furthermore, even if a large 
core is formed, a stellar density cusp can regrow via en¬ 
ergy exchanges between stars moving in the gravitational 
field of the MBH remnant. Black holes with masses of 
1 — 5% that of the NS C, such as those considered in 
IBekki fc CrahanTI ()201(lf) . correspond to low mass galax¬ 
ies with a short (< 10® yr) relaxation time. From equa¬ 
tion (IMl) one finds that such a low mass ratio corresponds 
to a MBH mass of ~ 5 x 10® Mq. Using equation (l3^ . 
this corresponds to a host stellar spheroid velocity dis¬ 
persion of Ri 50 km/s. In such low mass galaxies, the 
cusp will regenerate itself in about 1 Gyr from the infall 
of a second black hole, as collisional relaxation d rives the 
stella r distribution toward its steady state form (IMerritt! 
1201311 . We conclude that in low mass spheroids the effect 
of MBH binaries on the host NSC density distribution is 
expected to be minimal, and their imprint on the nuclear 
properties is expected to be erased after a short time. 

The probability that a galaxy of the size of the Milky 
Way experienced a major merger (i.e., with mass ratio 
large r than 0.3) after z = 2 is about 50% (|Fakhouri et al.l 
[2013). Galaxies of such luminosities have relaxation 
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times that are longer than the age of the galaxy, but 
which are short enough that two-body relaxation might 
have affected their central clusters. These galaxies might 
have nuclei that are still evolving toward their steady 
state after they were dynamically heated by a MBH bi¬ 
nary. Thus, the results of our models imply that “under- 
massive” NSCs in galaxies with bulge mass Mbuige ~ 
10 ^*^ Mq might still carry an imprint of the merger his¬ 
tory of their host galaxy, in the form of a density core that 
extends out to a small fraction of the central MBH influ¬ 
ence radius. Whether these NSCs will turn out to harbor 
parsec-scale cores will be addressed in the future with 
high resolution imaging, which will en able to resolve the 
crow ded stellar environment of NSCs (|Gullieuszik et al] 

[Ml . 


7. CONCLUSIONS 

NSCs and MBHs represent important central compo¬ 
nents of galaxies that bear witness to their formation 
history. Our aim has been to generate a numerical frame¬ 
work to better understand how the central structures of 
galaxies formed, and how the evolution of MBHs and 
NSCs is connected to that of their host galaxies. 

We have studied the evolution of NSCs in a cosmolog¬ 
ical context, by taking into account the growth of MBHs 
by merging subhalos containing both MBHs and NSCs. 
A semi-analytical galaxy formation model is applied to 
follow the evolution of dark matter halos along merger 
trees, as well as that of the baryonic components. Specifi¬ 
cally, the evolution of dark matter halos along the merger 
trees includes their baryonic components such as hot gas, 
stellar and gaseous bulges, and stellar and gaseous galac¬ 
tic disks. We study the simultaneous evolution of MBHs 
and NSCs in our galaxy formation model. The main re¬ 
sults of our work are summarized in the following: 

(1) The mass growth of NSCs is seen to be regulated 
by inspiral of star clusters, in-situ star formation as well 
as galaxy and MBH mergers. We have found that both 
cluster inspirals and in-situ star formation contribute a 
significant fraction of the total mass of NSCs. 

(2) We found that in-situ star formation (as well as 
growth through migrating clusters), generates NSC - host 
galaxy scaling relations that are much shallower than the 
same correlations for MBHs, in agreement with observa¬ 
tions. 

(3) In our model, the merger history of galaxies after 
z ~ 2 plays a key role in determining the structure of 
the NSCs observed today. Core-depletion due to grav¬ 
itational slingshot of host galaxy stars by inspiralling 
MBHs forming during galaxy mergers has a negligible 
impact on the evolution of NSCs in low mass galaxies, 
while it leads to their full disruption in galaxies more 
massive than ^ 10^^ Mq. In galaxies of intermediate 
luminosity 10^° Mq), MBH mergers cause the partial 
dissolution of the clusters. Such partially eroded clus¬ 
ters appear at z = 0 significantly underweight relative to 
linear NSC-host spheroid scaling correlations. The pres¬ 
ence of an under-massive NSC population results in an 
increased scatter of the NSC-host galaxy empirical cor¬ 
relations at high galaxy luminosities, for which we find 
evidence in observational data. 

(4) We identify Mmbh ~ 10®Mq as the critical value of 


MBH mass above which NSCs are significantly affected. 
Once the galactic MBH mass grows above this value, 
MBH binary mergers become highly disruptive for the 
central NSC leading to the formation of a central low 
density core. For galaxies more massive than the Milky 
Way, after such low density core forms, the central galac¬ 
tic regions are likely to remain depleted of stars through¬ 
out the subsequent evolution of the galaxy. 

(5) Our models predict that the fraction of nucleated 
early type galaxies containing an accreting MBH (with 
bolometric luminosity L > 10^° Lq) increases from 5% 
at Mgx ~ 10® Mq to 30% at Mgx ~ 10^^ Mq. These frac¬ 
tions are found to be much smaller for late type galax¬ 
ies, for which only ~ 10 % of nucleated galaxies with 
Mgx ~ 10^^ Mq also contain an active MBH. Among 
galaxies of all types, the fraction of mixed systems con¬ 
taining both a MBH and a NSC is about ^ 20% at 
10® Mq < Mgx < lO^i Mq. 

(6) We found that the results shown in this paper are 
quite robust when we change the cluster formation effi¬ 
ciency, /gc ^0.2; when we allow /gc to vary with galactic 
properties (e.g. we set it to 0.07, 0.04 and 0. 5 in disk, 
quies cent and starburst galaxies respectively; IKruiissei] 
l2012fl : when we change the MBH seed model (i.e. light- 
seed models vs heavy-seed models, with several halo oc¬ 
cupation numbers at high redshift) as well as other de¬ 
tails of our model, such as merger-tree resolution, initial 
redshift of the simulations, prescriptions for AGN feed¬ 
back, etc. 

In conclusion, our study supports a scenario in which 
black holes form inside NSCs with a low-mass fraction, so 
that NSCs are initially the dominant central component 
of galaxies. After their formation, NSCs and MBHs grow 
in parallel on their own scaling correlations with host 
galaxy properties, with NSCs dominating in low mass 
spheroids and MBHs dominating in high mass galaxies. 
The simple fact that NSCs and MBHs grow on different 
scaling correlations explains the well known transition 
from MBH- to NSC-dominated galaxies as one proceeds 
from dwarfs to giant ellipticals, without the need of in¬ 
voking competitive feedback processes from young NSCs 
and/or AGN activity. 
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